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FOREWORD 

In  the  present  summary  report  an  account  is  given  of  theories  and 
experiments  designed  to  investigate  the  irregular  structure  of  the 
ionosphere  by  studying  the  correlation  between  radio  waves  reflected 
at  closely  spaced  frequencies.  In  part  I  of  the  report  several 
different  models  are  constructed  of  the  ionosphere,  and  the  relation 
between  observable  quantities  on  the  ground  and  the  properties  of 
the  ionospheric  irregularities  are  derived.  In  part  II  of  the 
report  a  set  of  experiments  on  reflections  from  the  F-layer  at 
closely  spaced  frequencies  is  described  and  the  results  of  the  obser¬ 
vations  are  discussed  in  relation  to  the  results  of  part  I. 

Part  III  gives  detailed  description  of  the  equipment  used  in  the 
observations. 

The  work  described  is  sponsored  by  the  US  Air  Force  through  Office 
of  Aerospace  Research  under  contract  AF-61(052)-500, 
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Part  I 


ON  THE  DESCRIPTION  OF  ELECTROMAGNETIC  WAVES  REFLECTED  FROM  A  RANDOM 
MEDIUM 


SUMMARY 


Electromagnetic  waves  reflected  from  a  random  medium  are 
described  in  terms  of  spatial  correlation  and  of  correlation 
between  signals  reflected  at  adjacent  frequencies.  The  pro¬ 
pagation  of  these  correlation  functions  are  studied  both  for 
limited  and  unlimited  reflectors.  Several  different  reflec¬ 
tor  models  which  could  give  rise  to  the  random  properties  of 
the  reflected  waves  are  then  discussed  and  various  parameters 
are  related  to  the  observable  correlation  functions. 


1  INTRODUCTION 

Electromagnetic  waves  in  passing  through  or  being  reflected  from  a  spa¬ 
tially  random  medium  such  as  the  ionosphere  are  in  many  cases  perturbed 
in  a  more  or  less  random  manner.  A  single  plane  wave  passing  through 
such  a  medium  will  often  lose  its  plane  wave  character  and  develop  into 
an  angular  spectrum  of  waves.  Similarly  a  short  pulse  passing  through 
such  a  medium  will  be  widened  out  because  of  the  irregularities.  Often 
the  medium  might  also  vary  with  time  and  the  wave  passing  through  be¬ 
comes  time  varying.  This  means  that  a  wave  which  initially  consists 
of  a  single  spectral  line  gradually  will  develop  into  a  frequency  spec¬ 
trum.  All  of  these  perturbations  may  also  be  described  in  a  different 
way.  A  wave  which  is  initially  coherent  over  a  large  spatial  area 
will  lose  this  coherence  as  the  wave  is  progressing  through  a  random 
medium.  The  definite  phase  and  amplitude  relationship  initially  ex¬ 
isting  in  a  modulated  wave  will  gradually  vanish  with  the  wave  travel¬ 
ling  through  a  random  medium.  With  a  time  varying  medium,  the  rela¬ 
tion  of  fields  separated  by  a  time  interval  will  become  washed  out  as 
the  wave  penetrates  into  the  random  medium. 
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The  random  nature  of  the  diffracting  medium  will  only  permit  a  statis¬ 
tical  description  of  the  components  of  the  wave-fields.  The  different 
field  components  will  be  functions  of  position  in  space  r  ,  of  time 
t  and  of  the  frequency  u  of  the  applied  signals.  In  this  report 
we  describe  the  field  components  by  means  of  a  covariance  relating  the 
fields  observed  at  spatial  distances  A?  and  at  frequencies  separated 
by  Au.  In  what  follows  we  mostly  leave  out  the  description  of  the  ef¬ 
fects  of  the  random  time  variation  and  thus  only  deal  with  a  covari¬ 
ance  being  a  function  of  (4r,iw). 

For  signals  reflected  from  the  ionosphere  the  space  and  time  correla¬ 
tions  have  been  extensively  studied  (1,  2,  3).  Often  the  interpre¬ 
tation  is  rather  difficult  in  that  the  size,  shape  and  location  of  the 
irregularities  in  the  medium  giving  rise  to  the  random  field  cannot  be 
deduced  without  a  considerable  amount  of  additional  hypotheses.  For 
signals  transmitted  through  the  ionosphere  from  radio  stars  or  earth 
satellites  the  interpretation  of  the  wave  structure  is  considerably 
less  ambiguous  (1,  4,  5,  6).  Important  data  have  been  obtained  on 
ionospheric  irregular  structure  from  such  observations.  In  view  of  the 
many  unsettled  questions  concerning  the  reflection  mechanism  and  the 
irregular  structure  in  the  ionosphere  ,  it  appears  that  new  observa¬ 
tional  techniques  should  be  attempted  in  order  to  provide  additional 
or  supplementary  information.  One  of  these  techniques  might  be  the 
correlation  of  signals  reflected  at  spaced  frequencies,  or  alterna¬ 
tively,  the  study  of  the  distortion  of  extremely  short  pulses  reflec¬ 
ted  from  the  ionosphere.  Attempts  in  this  direction  have  been  made 
before,  but  the  interpretation  of  previous  data  might  be  of  doubtful 
validity  (7). 

In  this  part  of  the  report,  therefore,  we  want  to  study  theoretically 
the  relationship  between  space  and  frequency  correlation  functions 
for  several  different  models  of  the  irregular  structure  of  the  iono¬ 
sphere.  We  are  in  particular  attempting  to  examine  to  what  extent  the 
frequency  correlation  technique  could  provide  additional  information  on 
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the  ionospheric  structure  and  under  which  circumstances  this  be  the 
case.  Section  2  below  therefore  deals  with  some  general  properties  of 
the  statistical  functions  of  the  wave-field.  Section  3  discusses  the 
relation  between  the  properties  of  the  wave-field  and  the  structure  of 
the  irregular  medium  for  a  number  of  different  models. 


2  DESCRIPTION  OF  THE  DIFFRACTED  FIELD 


In  this  section  we  give  a  brief  description  of  the  properties  of  the 
random  diffracted  field.  We  assume  that  the  field  is  given  over  a 
screen  at  z  =  o  and  that  field  variations  only  take  place  in  the 
x-direction,  hence  the  field  is  one-dimensional.  We  furthermore  as¬ 
sume  that  a  statistical  ensemble  of  such  random  screens  are  given  and 
that  the  statistical  properties  of  the  ensemble  is  given.  The  actual 
physical  mechanism  giving  rise  to  the  random  field  distributions  over 
the  screen  is  not  discussed  here.  Certain  simple  models  are  dealt 
with  in  the  next  section.  Knowing  the  statistical  properties  of  the 
field  in  the  plane  z  -  0  we  inquire  about  the  propagation  of  these 
properties  into  the  half  plane  z  >  0.  This  approach  has  been  used  by 
J  A  Ratcliffe  (8)  and  it  has  been  extended  by  others.  Our  main  excuse 
for  repeating  some  of  the  already  known  arguments  here  is  that  we  also 
want  to  include  in  our  discussion  the  relation  of  two  harmonically 
varying  fields  at  slightly  different  frequencies.  The  two  cases  of  an 
infinitely  extended  screen  and  that  of  a  spatially  limited  one  are  con¬ 
sidered  separately. 


2.1  Covariance  for  infinite  screen 


Let  the  field  over  the  screen  at  z  -  0  and  at  a  frequency  u  be  de¬ 
noted  by  g(x,0,w).  This  field  can  be  expressed  in  terms  of  an  angu¬ 
lar  spectrum  F(S,w)  through  (8): 


'■a 
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corresponds  more  closely  to  the  actual  physical  problem  under  study, 
see  also  (9). 

The  nature  of  the  random  field  will  now  be  described  in  terms  of  the 
complex  correlation  coefficient  R(Ax,Aw,z)  defined  by: 

R(AX,pu,z)  =  ^g*(x,z,u)  g(x+4X,z,iu+Au)^av  (2.4) 

where  the  asterisk  denotes  the  complex  conjugate.  At.  the  screen  the 
statistical  properties  should  be  obtained  by  putting  z  =  0. 

In  order  to  relate  the  correlation  at  the  screen  to  the  correlation  at 
some  point  in  front  of  the  screen  we  observe  that: 


g(x,z,«)  =  j  dS  |  d£  e({,0,u)  exp{  i  £  (S(x-£)+Cz)J  (2.5) 

which  follows  from  a  combination  of  formulae  (2.1)  -  (2.3).  We  there¬ 
fore  find: 


R(Ax,pu,z)  =  iU2  ffff  dSdS'd{dt'(g*(5»0>ui>g(? '  ,0,u+au)\  . 

(2„c)  1111  \  /av 

•  exp{  -i  £  (S(x-$)+Cz)  +  i  (S' (x+ax-{ '  )+C' z) ) 

(2.6) 

Because  the  statistical  properties  of  the  field  at  the  screen  are  in¬ 
dependent  of  position,  the  average  under  the  integral  sign  only  depends 
on  the  difference  p  *  {'  -  5.  Integration  over  6  therefore  gives: 

R( Ax, Au, z )  =  w(w*AM^  f  f  [  dSdS'dp  R(  p ,  Aw,  0 ) .  A  (—  -  SlAtt  S')  • 

2it*o2  c  c 

expf  -i  “  (Sx+Cz)  +  i  “’^“(S'lx+Ax-e)  +  C'z)) 


(2.7) 
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where  i(  )  is  a  Dirac  delta-function.  Different  components  in  the 
angular  spectrum  therefore  contribute  at  the  two  frequencies,  viz.: 


(2.8) 


Carrying  out  the  integration  over  S'  in  (2.7)  leaves  us  with  a  double 
integral  for  the  correlation  coefficient: 


R(Ax,Aw,z)  =  - ' —  [  dS  exp{  i—  S’Ax  ♦  i"r(C' (w+Aw)  -  C'w)  } 

2*c  1  c  u 

•  |  dp  R(p,Aw,0)  exp(-iSp  ®)  (2.9) 

Let  us  now  examine  two  special  cases,  the  first  being  Aw  =  0.  In  this 
case : 


R(  Ax ,  z )  =  II  dSdp  R(  p  ,  0 , 0 )  exp(  i|  S(Ax-p)}  =  R(Ax,0,0) 

(2.10) 

The  spatial  correlation  therefore  remains  invariant  with  increasing 
distance  from  the  screen,  as  is  well  known  (8).  In  the  second  case 
Ax  =  0  and  we  obtain: 


R(Aw,z)  =  II  dSdp  R(p,Aw,0)  exp{-i^  Sp  +  i~(C  (w+Aw)-Cw)  J 

(2.11) 

It  can  be  seen  that  this  expression  does  vary  as  the  distance  z  is 
increased.  The  explanation  for  the  variation  of  frequency  correlation 
with  distance  z  is  best' seen  as  follows.  With  increasing  z  the 
exponential  will  become  more  and  more  rapidly  oscillating  function  of 
Au.  Eventually  contributions  will  only  arise  when  Aw  :  0  in  the 
correlation  function,  hence: 
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R  ( Am ,  z  >  -*■  J  d  S  exp{i~(C'  (ui+Aw)-Cw) }  •  j  dp  R(p,0,0)  expl-i^fsp} 

(2.12) 

The  last  part'  of  this  is  equal  to  the  angular  power  spectrum  at  the 
frequency  w.  Expanding  furthermore  the  first  exponent  to  first  order 
in  Aw  we  find: 


R(Au,z)  -  f  dS  /|F(S,w)|2)  exp{  i  — =j=T  }  (2.13) 

)  \  '  av  /,  .2 

c  y  i  -  s 

In  this  asymtotic  form  the  result  is  fairly  easy  to  understand  if  it 
is  rederived  in  a  different  way  by  making  use  of  the  relation  between 
the  average  power  pulse  response  (or  "delay  spectrum")  D(t)  and  the 
correlation  function  R(aw).  Here  t  represents  the  delay  of  the 
network  which  is  equivalent  to  the  propagation  circuit.  The  relation 
is  (10): 


R(Aw)  =  |  D(r)  exp{  iAw*T  }  di  (2.14) 

which  is  valid  as  long  as  (Aw/u)  <<  1.  This  is  also  the  condition 
for  (2.13)  to  be  valid. 

In  Figure  2.1  it  is  seen  that  the  delay  t  is  given  by 


Substitution  into  (2.14)  gives: 


R( Au | z ) 


dSif 


yr 


D(- 


v/T 


-)}  exp(i 


(2.16) 
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Identification  of  the  expression  within  the  curly  brackets  with  the 
angular  power  spectrum  brings  this  into  agreement  with  (2.13).  The 
interpretation  is  now  clear:  with  increasing  distance  from  the  screen 
and  with  an  opening  angle  in  the  power  spectrum  which  does  not  change 
with  distance  the  spread  in  time  delays  is  increasing  and  the  frequency 
correlation  range  narrows.  At  the  same  time  this  suggests  that  with  a 
source  which  is  limited  in  extent  the  frequency  correlation  might  not 
change  much.  This  is  examined  in  detail  below. 


2 . 2  Covariance  for  finite  screen 


Let  us  next  assume  that  the  random  diffraction  fields  are  viewed 

.  .  W  W 

through  a  slit  at  z  =  0  which  extends  from  x  =  -  t;  to  x  =  y  ,  and 

that  the  statistical  properties  are  independent  of  position  within  the 

slit.  We  furthermore  assume  that  the  scale  of  the  irregularities  is 

much  smaller  than  the  width  of  the  slit.  Under  these  assumptions  we 

find  that  equation  (2.6)  must  be  replaced  by: 


2 

RtAXjAu.z)  =  j  2  |  ^d(JdpR(p  ,Au,0)  JJ  dSdS'  exp( -i^(  S  (x-C )  +Cz )  + 

"7 

+  j2i^(S'(x+6x-6-p)  +  C '  z ) )  (2.17) 

Because  of  the  limited  extent  of  the  source  the  fields  will  die  out 
with  distance  from  the  screen  and  the  above  correlation  function  will 
not  be  normalized  as  in  the  case  of  an  infinite  screen.  Integration 
over  S  and  S'  gives  the  unnormalized  correlation: 

1  /—  —■  —  .  f 

R(Ax,Aw,z)  s  TTc^z  vw(w+Aw)  expji— 1 “ — }  Jdp  R(p,Aw,0) 

exp{i2cz(ix"B>2)  I  d5  expti«C  <  Ax-P  >+-|^Cx-e  )2 } } 


(2.18) 
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The  interpretation  of  this  is  most  convenient  if  we  consider  the  same 
two  special  cases  as  above,  viz,  those  of  ax  -  0  and  Au  -  0.  When 
Au  =  0  : 


R(  Ax , z ) 


(2.19) 


Here  we  see  that  when  w  ■*  «  with  z  constant  the  sine  term  will  tend 
to  a  4-function  and  contributions  to  the  last  integral  will  only  arise 
where  p  =  Ax  ,  and  it  will  be  found  that  R(Ax,z)  will  be  constant 
and  the  same  as  at  the  screen.  When  we  let  z  ■*  »  with  W  constant, 
the  integral  will  become  a  more  and  more  slowly  varying  function  of 
Ax.  The  distance  in  space  over  which  the  correlation  remains  high 
therefore  increases  with  distance  from  the  screen.  Note  that  this 
applies  to  a  finite  screen  and  with  ensemble  averaging! 


When  Ax  *  0: 


R(Aw,z)  =  v4(w+Aw)  exp(  i  A“c  Z  (1  +  i(^)2)}  • 

•  |  dp  R(p,Aw,0)  exp(  i  (p2  -  2xp ) }  (2.20) 

In  performing  the  integration  we  had  to  assume  that  W2  is  small  com- 

c 

pared  with  8z(~).  Here,  as  z  -►  »  ,  apart  from  the  rapidly  oscilla¬ 
ting  phase  factor  outside  the  integral  the  frequency  correlation  tends 
to  become  constant  and  independent  of  z. 

The  conclusion  from  this  section  can  therefore  be  stated  briefly  as 
follows : 

For  an  infinite  screen  the  space  correlation  is  independent  of  distance 
from  the  screen  whereas  the  frequency  correlation  narrows.  For  a  fini¬ 
te  screen  the  space  correlation  widens  out  with  distance  from  the 
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screen  whereas  the  frequency  correlation  tends  towards  an  asymptotic 
value  with  z  only  determined  by  the  properties  of  the  correlation  cn 
the  screen. 

An  extension  to  a  two-dimensional  screen  is  simple  in  principle,  but 
it  is  not  carried  out  here  because  it  does  not  appear  that  any  essen¬ 
tially  new  points  of  physical  interest  are  brought  out  thereby. 


3  THE  ORIGIN  OF  THE  RANDOM  WAVE  FIELD 


In  the  previous  section  we  assumed  the  field  variations  to  be  described 
statistically  at  some  plane  z  =  0.  We  then  went  on  to  study  the  pro¬ 
pagation  of  these  properties  away  from  the  screen.  -  In  the  present 
section  we  must  inquire  about  the  different  physical  mechanisms  which 
could  give  rise  to  the  random  field  variations  over  the  reference  screen. 
Because  we  want  to  consider  waves  reflected  from  the  ionosphere  we  do 
not  discuss  the  problem  of  a  wave  being  transmitted  through  an  irregu¬ 
lar  medium  which  would  be  appropriate  in  the  case  of  radio  star  or  earth 
satellite  studies.  The  problem  of  reflection  from  an  irregular  medium 
is  so  complex  that  we  must  try  to  construct  models  which  are  simple 
enough  to  be  discussed  without  too  much  difficulty,  yet  sufficiently 
detailed  to  convey  some  of  the  features  of  the  actual  reflection  mecha¬ 
nism.  Two  types  of  models  will  be  considered.  The  first  type  is  a 
surface  model.  The  waves  are  regarded  as  being  reflected  from  randomly 
corrugated  surfaces  with  reflection  coefficient  unity.  The  second  type 
is  a  volume  mechanism  where  irregularities  in  electron  density  are  em¬ 
bedded  in  the  ionospheric  layer  below  the  reflection  level.  The  irre¬ 
gularities  are  considered  as  weakly  scattering. 

For  simplicity  we  consider  only  an  incident  plane  wave  at  vertical  in¬ 
cidence  on  the  ionosphere.  An  arbitrary  incident  wave  can  be  construc¬ 
ted  from  an  angular  spectrum  of  such  plane  waves. 
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3.1  Reflection  from  a  single  corrugated  surface 


Figure  3.1  Reflection  from  corrugated  surface 


Consider  the  surface  shown  in  Figure  3.1.  The  random  corrugations  are 
in  the  x-direction  only  and  the  function  h(x)  describes  the  devia¬ 
tion  of  the  surface  from  its  mean  position. 

Provided  the  surface  gradients  are  not  too  steep  and  the  structure  of 
the  surface  too  small,  the  field  referred  to  z  =  0  which  is  the  mean 
position  of  the  screen  becomes: 

g(x ,0 ,u)  =  exp{  -  i  “  •  2h(x)  }  (3.1) 

It  is  implied  that  the  incoming  plane  waves  are  of  unit  amplitude  and 
that  the  reflection  coefficient  of  the  surface  is  unity.  The  corre¬ 
lation  function  at  the  screen  becomes: 

R(Ax,Aw,0)  =  ^exp{  i  £(u>h(x)  -  (w+Aw)  h(x+Ax))}^av  (3.2) 

In  order  to  make  some  progress  with  this  it  will  be  assumed  that  h(x) 
and  h(x+Ax)  are  distributed  in  accordance  with  a  joint  normal  dis¬ 
tribution  with  the  following  parameters: 
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(h<x)2)av  =  (h(x+4x)2)av  =  hm2 

^h(x)^)av  =  ^h  (x+Ax)^av  =  0  (3.3) 

^h(x)h(x+ix)^  av  =  hm2  p  ( Ax ) 

With  these  assumptions  it  follows  that: 


R(Ax,Aw,0)  =  exp(- (~~)  2  ( w(  w  +  Aw)  ( 1  -  p  ( Ax ) )  +  j  )}  (3.4) 

For  Aw  =  0  this  simplifies  to 

2whm 

R(Ax,0,0)  =  exp(  -  ( — ”  )  (1  -  p ( Ax ) ) }  (3.5) 

and  this  is  the  same  as  the  correlation  at  the  ground  in  the  case  of 
ionospheric  reflections  provided  the  aerial  polar  diagram  is  not  so 
narrow  that  the  screen  appears  as  if  limited.  The  result  (3.5)  is 
well  known  from  previous  investigations  (11).  We  now  have  to  decide 
upon  a  specific  form  of  the  correlation  function  p(Ax)  at  the  re¬ 
flecting  surface  before  further  progress  can  be  made.  For  lack  of 
specific  information  we  make  the  convenient  choice 

p ( Ax )  =  exp(-Ax2/l2)  where  1  is  the  scale  of  the  corrugations.  The 
joint  correlation  function  at  the  receiver  then  becomes,  see  (2.9): 


R(Ax,Aw,z)  *  exp(-  2 ( — 5 — )z)  J  dS 

exp{i  *  S • Ax  +  i  §(C'(w+Aw)  -  Cw)}  • 

•  |  d?  exp{-(— £^) 2 ( w( w+Aw) (1  -  exp(-£2/l2) ) ) )  • 

exp(-  i  S?  £)  (3.6) 


We  now  distinguish  between  the  case  of  a  deep  reflector,  or  hm  >  A 
and  a  shallow  reflector,  hm  <  A . 
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It  is  convenient  here  and  in  what  follows  to  introduce  the  r.m.s. 
opening  angle  8m  of  the  angular  power  spectrum  to  simplify  the  for¬ 
mulae.  In  the  case  of  a  deep  reflector  we  find  approximately  that: 


8  <-fr 


(3.7) 


Introduction  of  this  in  formula  (3.6)  gives,  still  for  a  deep  reflec¬ 
tor: 


R(Ax,Aw,z)  = 


exp 


{  i*|S  - 


J\  -  i  8  2*law(l+^S) 

v  wi  r>  i.k 


expt 


6  2  Ax2 


1  -  i  8  2  -4u(l*— ) 
me  u 


(3.8) 


When  flu  *  0  this  reduces  to: 


R(Ax,z)  1  exp(  -  i  e^2  (^)2  Ax2 


(3.9) 


which  is  also  a  reasonable  approximation  for  R(Ax,0)  in  this  parti¬ 
cular  case. 

When  Ax  =  0  and  when  phase  factors  are  neglected  we  have: 


R(Aui,z) 


/ 


1  -  i  8  2  iAwd+iSL) 
me  w 


exPl  -  hm2 


(3.10) 


The  first  factor  which  is  due  to  the  angular  spectrum  of  waves  com¬ 
prising  the  total  signal  drops  to  half  of  its  magnitude  at  Au  =  0 
when: 


8m2  f  ».  -  /r 
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or  when: 


Aw 


z 


(3.11) 


The  exponential  factor,  on  the  other  hand  drops  to  about  one  half  when 


2(Aii)2.h  2  =  i 

c  m 

or  when 

Awc  =  (3*12) 

c  -n.  % 


For  this  model  it  is  well  known  that  the  observation  of  spatial  corre¬ 
lation  only  gives  information  on  the  ratio  (hm/l).  If  the  frequency 
correlation  is  to  provide  additional  information  we  must  require  that: 


Aw q  <  Awz 


which  means  that: 


In  the  other  case,  that  of  a  shallow  reflector,  the  r.m.s.  opening 
angle  for  the  random  part  of  the  spectrum  becomes: 


(3.14) 


The  joint  correlation  function  at  distance  z  from  the  screen  in  this 
approximation  again  from  formula  (3.6)  becomes: 
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R(AX,Aw,  z) 


exp( 


AuZ 


-  2  (- 


A  wh 


12.2 


(1 


2h„  5 

(”c  )  w(u>  +  Aui)  { 1  - 


expl-ty1)2  -  — ■1-— ~  ) 

1-10  2  rAw 
m 

A  .  i  i  ! 

v  me 


-}) 


(3.15) 


which  for  Aw  =  0  reduces  to: 

2hmw  aj. 

R(Ax,z)  *  (l  -  ( — " — )  (1  -  exp(  -  (— )2  })}  (3.16) 

We  therefore  see  that  the  correlation  does  not  drop  to  zero,  but  to  a 
value  determined  by  the  depth  hm  ,  and  that  the  width  is  determined 
only  by  the  scale  of  the  corrugations  1.  When  Ax  =  0  we  obtain: 

A  whm  A  w  , 

R(Aw,z)  *  exp{  -  2  (  “ — )  *  i  ~  z  ] 

2hm  , 

{l  -  (-^)  w(w+Aw)fl -  *-  ■—--}]  (3.17) 

J \  -  i  0  2  ~  Aw 
v  me 

Aw  .  uhm 

Because  —  <<  1  and  since  ~~£~  <  1  by  assumption,  the  first  expo¬ 
nential  factor  is  very  close  to  unity.  The  width  of  the  random  part 
of  the  returned  waves  is  simply  determined  by  the  width  of  the  angu¬ 
lar  power  spectrum,  and  no  more  information  can  be  obtained  from  fre¬ 
quency  correlation  than  from  space  correlation. 


3 . 2  Reflection  from  two  corrugated  screens 

We  next  imagine  the  reflection  at  the  two  frequencies  to  take  place  at 
different  surfaces  corresponding  to  the  deeper  penetration  into  the 
ionosphere  by  the  higher  frequency  wave.  We  make  the  very  crude  as¬ 
sumption  that  both  waves  travel  undisturbed  up  to  their  reflection 
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level,  and  that  the  reflection  occurs  with  a  reflection  coefficient  of 
unity.  We  furthermore  take  the  electron  density  in  the  mean  to  be  a 
linear  function  of  height  z  ,  viz.: 


N ( z )  =  -  g(z  -  zQ)  z  <  z0 

g  >  0 


(3.18) 


(note  that  z  is  positive  downwards!).  The  difference  in  mean  height 
of  reflection  caused  by  the  frequency  change  Aw  becomes: 


2Me0 

Az  *  j  w  •  Aw  =  -  A  Aw  (3.19) 

ge 


Each  of  the  two  reflecting  surfaces  are  regarded  as  being  described 
statistically  in  the  same  way  as  the  single  surface  in  the  previous 
section,  and  both  with  the  same  statistical  parameters.  The  two  sur¬ 
faces  are  not  identical,  but  differ  more  and  more  with  increasing  mean 
separation  Az.  In  this  case  we  let  the  deviations  from  the  mean 
heights  be  denoted  by: 


h(x,z)  at  height  z 

h(x,z+Az)  at  height  z*Az 


(3.20) 


and  we  take  the  relation  between  the  height  deviations  to  be: 


\ 


(3.21) 


The  situation  is  illustrated  in  Figure  3.2  showing  three  different 
surfaces.  For  convenience  the  reference  plane  is  now  chosen  at  the 
mean  height  of  reflection  for  the  lower  frequency. 


For  plane  waves  of  unit  amplitude  we  obtain  for  the  reflected  fields 
at  the  reference  plane  z  =  0, 
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h(x,  az2) 


Figure  3.2  Reflection  from  two  corrugated  surfaces 


g(x,0,iu)  =  expj  -  i  ^  2  h(x,0)  f 


end 


w+Au)  , 

g(x,0,(u+4u)  =  exp(  -  i  — 2Ch(x,Az)  -  a z)  }  = 

.  ai  +  A  b) 

=  exp{  -  i  — " —  2(h(x,Az)  *  A-Au))  (3.22) 

For  the  complex  correlation  at  the  reference  plane  we  now  obtain: 

(u+Au  , 

R(Ax,Aw,0)  =  exp{  -  2i  — " —  A z  J  • 


exp{  -t“~)  2  ( w(w+Aw)  (1-p  (Ax.AAw) )  +  (3.23) 

Again  we  have  to  make  an  assumption  about  the  correlation  function 
p(Ax,Az)  in  order  to  make  further  progress.  We  arbitrarily  put: 


exp{  -  ( -  + 

1» 


P (AX,Az) 


(3.24) 
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where  m  must  be  interpreted  as  a  vertical  scale.  Distinguishing 
again  between  deep  and  shallow  reflectors  as  before  one  finds: 

For  deep  reflector: 


,  .  i  .  m±Am  Am* z. . 

R(Ax,Aw,z)  =  exp(  i  (2  — ~  Az  ♦  — ” )} 


expl  -  2>)2(1  ♦  2  &***!,  iu2  | 


/ 


1  -  i  9  2  i  Au(l  +  A“) 
me  w 


exp{ 


$  x2  .  e  2  .  w(w+Aw) 

_ m  c  2 

1  -  i*m**fAu(l+~) 


(3.25) 


When  Aw  s  0  we  of  course  must  obtain  the  same  as  in  the  previous 
section  because  the  two  reflection  surfaces  then  coincide.  When 
Ax  =  0  we  have,  neglecting  phase  factors: 


R(Au,z)  * 


-  i9  2  £aui(i  +  Ail) 


i  -  A2w(m+Aw)  . 

exp{  -  2 ( — )'(1  +  2  -  )  Aw2  } 

*  C  *»  ’ 


(3.26) 


To  interpret  this  we  note  that: 


2 


Mco 


w 


L 


(3.27) 


where  we  have  introduced  the  plasma  frequency 
defined  by: 


P 


and  the  length  L 
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L 


-1 


1 


N 


dN 


dz 


(3.28) 


At  the  reflection  level  the  frequency  equals  the  plasma  frequency  wp 
and  A  therefore  simplifies  to: 


Substitution  into  the  exponential  term  in  (3.26)  gives: 


Because  L  is  the  half  thickness  of  the  layer  and  because  m  is  the 

L 

vertical  scale  of  the  irregularities  we  expect  ~  >>  1.  The  half  width 
of  the  exponential  factor  now  becomes: 


Au 


c 


1 

4 


c 

hT 


(3.29) 


which  must  be  compared  with  the  width  of  the  angular  spectrum  factor 
which  is  the  same  as  for  a  single  reflector,  see  (3.11).  The  condition 
for  the  vertical  scale  to  be  detectable  becomes  4uc  <  A«z  which  leads 
to: 


2  ‘  4  •  S*  '  m  *  7*2  =  <3.30) 

m  m 

This  result  should  be  compared  with  (3.13).  Since  8  is  of  the  order 
10-1  and  since  the  half  thickness  of  the  F  layer,  say,  is  about  half 
the  reflection  height  it  seems  quite  likely  that  on  this  model  the  fre¬ 
quency  correlation  technique  could  provide  information  in  addition  to 

the  one  provided  by  spatial  correlation,  viz.  on  (“)  . 

"m 

Turning  now  to  the  case  of  shallow  reflectors  we  obtain  with  the  same 
approximations  as  in  the  previous  section: 
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f  w  +  AW  flllCS.  III.  2  l 

R(Ax,Aw,z)  :  exp{  i  <2  — - —  A z  +  — £ — )  -  2( — -  )  } 


All)'  z 


Awh, 


m  2 


2h™ , 

•(l-t-j2)  «(i»+4«>  {1  - 


Kpt-(f> 


1-iO  2  -  A«  m2 

- m__c - })  (3.31) 


Z1  -  i6m2  t  “ 


With  ax  s  0  and  with  phase  factors  suppressed  we  obtain: 


Awh. 


R  ( A  w ,  z  )  *  exp{  -  2 ( — T^)2  } 


.  A2Aw2 

2h  eXP("  2  > 

(1  -  (-SS)2u(u+A«.)  (1  -  . ZZ - 1 )  (3.32) 


/l  _  i6  2  -  Au 
v/  mo 


The  first  exponential  factor  is  essentially  unity,  and  it  is  the  second 
part  of  the  expression  within  the  last  curly  bracket  that  decides  the 
loss  of  correlation.  The  depth  to  which  the  correlation  falls  is.  again 
determined  by  the  depth  hm  of  the  reflector.  The  width  is  again  de¬ 
termined  by  two  competing  factors.  The  half  widths  of  these  are: 


7  ‘  t 


and 


Awz 


c 

z 


The  condition  for  the  vertical  scale  m  to  be  detectable  hence  becomes 


(3.33) 


where  *  is  the  free  space  wavelength.  If  om  is  of  the  order  of 
Yj  and  the  wavelength  is  not  smaller  than  the  vertical  scale  by  more 
than  a  factor  of  about  100  the  vertical  scale  might  be  detectable. 
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It  should  be  pointed  out  here  that  the  models  employed  so  far  are  of 
somewhat  doubtful  validity  in  that  the  waves  are  considered  to  be  tra¬ 
velling  as  if  in  free  space  up  to  the  reflection  level.  In  fact,  the 
wavelengths  are  increasing  toward  infinity  near  the  reflection  level 
and  it  might  be  that  the  vertical  scale  effect  on  the  frequency  corre¬ 
lation  is  largely  overestimated  by  this  method. 


•  3  Reflection  from  a  plane  reflector  in  a  blobbv  medium 

In  this  and  the  next  subsection  we  consider  the  random  amplitude  and 
phase  variations  as  imposed  on  the  waves  through  volume  effects  and 
not  as  above  through  the  reflection  from  a  rough  surface.  In  the  pre¬ 
sent  subsection  we  consider  the  reflecting  surface  to  be  plane  and 
perfectly  reflecting  and  to  be  embedded  in  a  blobby  medium.  The  ver¬ 
tical  distance  between  the  reflection  level  and  the  reference  level  at 
z=0  is  0  ,  and  all  the  scattering  blobs  are  considered  to  be  loca¬ 
ted  between  these  two  levels.  We  shall  limit  ourselves  to  large  scale 
blobs  by  assuming  the  horizontal  scale  to  satisfy  the  inequality: 


blobs 

2=0 


reflector 


°  O  O  o  0 
D  O  °  r 

°  o  o  <?C 


z 


Figure  3 


Reflector  in  blobbv  medium 
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i  >  y  d  •  x 

where  A  is  the  wavelength.  With  this  condition  satisfied  diffraction 
effects  within  the  blobby  medium  need  not  be  considered  and  the  propa¬ 
gation  problem  becomes  one  of  geometrical  optics.  We  furthermore  con¬ 
sider  the  irregularities  to  be  so  weak  that  partial  reflection  and 
bending  of  the  rays  can  also  be  neglected.  The  properties  of  the  medi¬ 
um  are  taken  to  vary  in  the  x  and  z-directions  only,  and  the  medium 
is  irradiated  by  plane  waves  along  the  negative  z-axis,  as  before. 

When  these  requirements  are  met,  the  random  field  over  the  reference 
plane  at  z  =  0  becomes: 


g(x,w,0)  =  exp{  -  2  i  j  K(x,iu,n)  dn  }  (3.34) 

o 

where  K(x,u,z)  is  the  complex  propagation  constant.  The  real  part 
is  related  to  phase  delay  and  the  imaginary  part  to  attenuation.  For 
the  unnormalized  correlation  function  at  the  reference  plane  we  ob¬ 
tain: 


R'  (Ax,Aw,0)  » 


^exp{2i  j  K  <x,w,n)dn  -J  K(x+Ax,w+Aw,n )dn)  . 

/  ,D(*i)  a  .D(w)  v 

5  ^exp{  2i(  j  Z)K(x,u,n)dn  -  tic-  J  K(x,w,n  )dn)  . 


(3.35) 


where  we  have  expanded  the  exponent  to  first  order  in  Au  and  where 
we  have  put  AK  =  K*(x,w,z)  -  K(x+Ax,u,z).  It  is  convenient  to  sepa¬ 
rate  phase  and  attenuation  effects  by  putting:  K  =  k  -  ix  ,  which 


means  that 
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AK  =  (k(x,w,z)  -  k(x+Ax,(u,z) )  -  i(x(x,w,z)  +  x (x*Ax,w ,Z ) )  = 

=  (k  -  k>)  -  i(x  ♦  x' )  (3.36) 

The  integral  J  (k  -  k')  dn  is  simply  the  difference  m  phase  delay 
1  o 

between  two  paths  separated  by  a  distance  A  x.  The  derivative  with 
respect  to  frequency  of  the  phase  delay  is  the  group  travel  time  T 
between  the  levels  z  =  0  and  z  =  -D.  In  the  absence  of  absorption 
the  following  simple  expression  therefore  applies: 

R(Ax,Aw,0)  =  ^xp{  2i  J  (k  -  k' )dn  -  iAu*T  ^av  (3.37) 

To  evaluate  the  space  and  frequency  correlations  at  the  reference 
plane  we  have  to  evaluate  the  probability  distribution  of: 


a  =  2  J  (k  -  k’)dn 


and 


8 


2(^  jVdn  -  £  dn)£ 


(3.38) 


The  neglect  of  absorption  effects  is  a  good  approximation  above  about 
100  km  for  frequencies  higher  than  about  one  MHz.  For  k  we  use  the 
expression: 


•  f  /TTY 


where  X  =  N^x,z£e 
Me0w 

sist  of  two  parts: 


(3.39) 


Here  the  electron  density  N(x,z)  con- 


N(x,z)  *  N0(z)  +  AN(x,z) 


(3.40) 


where  AN(x,z)  is  a  random  deviation  from  the  local  mean  value.  For 
the  random  part  we  make  the  assumptions: 
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(AN<*.*>)av  = 

(4N(x'z)2)av 

/aN(x+4x,z+Az) 


0 

N  2(z) 
m 

4N(x,z))av 


(3.41) 


N  2(z) 
m 


4Pi-(f)J 


- 


Because  a  and  8  can  be  considered  as  a  sum  of  a  fairly  large  num¬ 
ber  of  independent  variables  they  are  both  gaussian  with  mean  zero  and 
with  variances: 


2  /*"  <“)2  m(l-exp(-(*-)2})  f 


,D  Nm2 ( z ) 


dz 


(3.42) 


and 


yF  (fi)*  m 


,  -D  N  2  ( z ) 
i  f  _m _ 

2  I  6 

w  *  O  M 


dz 


(3.43) 


e  2 

where  B  s  Me  . ~  .  Again  we  distinguish  between  deep  and  shallow 
o 

phase  modulation  at  the  reference  plane. 


For  deep  modulation  the  spatial  correlation  becomes: 


R(Ax,0) 


*  exp(  -  \A T  (~)2  m  •  ♦  •  4x2  } 


(3.44) 


where 


♦  = 


fD  Nm2(z) 


dz 


(3.45) 


This  means  that  the  mean  square  angular  spread  becomes: 


yr  c^)2  cf )2  •  »  • 


(3.46) 


For  the  frequency  correlation  at  z  =  0 
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R( Au ,  0 )  :  exp{  -  J*  (^)2  m  •  »  •  (^)2  1 


(3.47) 


where : 


,D  N  2  (z) 


dz 


(3.48) 


The  half  width  of  the  corresponding  part  of  the  frequency  correlation 
becomes : 


Aid 


(3.49) 


whereas  the  half  width  caused  by  the  angular  spread  is  the  same  as  in 
(3.11).  The  condition  for  the  vertical  scale  m  to  show  up  therefore 
requires  that: 


z 


(3.50) 


From  this  it  appears  that  the  horizontal  scale  1  must  be  large,  the 
opening  angle  0m  small  and  the  ratio  7/*  large  for  the  vertical 
scale  m  to  be  observable.  If  the  irregularities  are  close  to  the 
reflection  level,  might  become  large  and  make  the  observation  of 

m  possible.  However,  the  model  does  not  appear  to  be  quite  satisfac¬ 
tory  near  the  reflection  level. 


For  shallow  phase  modulation  the  angular  width  of  the  random  part  of 
the  waves  becomes: 


(3.51) 


and  the  condition  for  the  non-angular  spectrum  part  to  be  the  dominant 
cause  of  the  loss  of  frequency  correlation  becomes: 
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w 


C 


v/m-t 


V57=B- 


<z) 


utn  \/m,T 


(3.52) 


Because  this  inequality  depends  very  strongly  on  the  irregularities 
near  the  reflection  level  through  t  we  shall  abstain  from  any  guesses 
as  to  the  validity  of  the  above  inequality. 


.4  Reflection  from  a  linear  gradient  in  electron  density  with  irregulari¬ 
ties  embedded 


The  scalar  wave  equation  has  recently  been  solved  by  a  perturbation 
analysis  for  weak  irregularities  embedded  in  a  linear  electron  density 
profile  (12,  13).  The  perturbation  analysis  only  holds  as  long  as  the 
perturbation  is  small  in  comparison  with  the  unperturbed  wave.  Put¬ 
ting: 


g(x,w,z)  =  go(x,w,z)  +  g,(x,w,z) 


(3.53) 


it  is  therefore  required  that  (g,/g0>  <<  1  .  This  limits  the  applica¬ 
tion  to  cases  with  shallow  fading,  or  to  cases  similar  to  those  of 
shallow  phase  modulation  at  the  reference  level  as  considered  above. 

As  we  have  seen  there  is  a  specular  component  present  in  this  case. 

The  joint  correlation  at  the  observer  becomes: 


^g*(x,<i,z)g(x+4x,w  +  flu,z)^  av  =  ^g0*(x,w,z)g0(x+Ax,w+Aw,z)^a, 

♦  ^g1*(x,u,z)»g, (x+Ax.w+Aw.z)^ 


(3.54) 


because  the  average  perturbation  is  zero.  The  first  term  does  not 
change  in  magnitude  with  Aw  and  Ax  ,  but  the  second  one  does.  It 
is  with  the  last  term  we  shall  be  concerned  here. 
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We  briefly  sketch  the  evaluation  of  g)(x,u>z).  The  scalar  wave 
equation  is: 


2  N  2 

{v2  +  k  2(1  +  — — 7  •  -r  •  z)  -  kQ2  — - — r  AN(x,z)}  g(x,z)  =  0 
Me0w  Me0w 

(3.55) 


this  equation 

all  quantities 

are 

expanded  into  Fourier  integrals  in 

x-direction! 

6l<x*z)  =  I 

dS  •  g((S,z)  • 

exp| 

[  i  S  kD  x  } 

(3.56) 

AN(x,z)  =  | 

dS  •  AN(S.z)  • 

exp| 

[  i  S  kQ  x  } 

...  .  2* 
with  kQ  =  y  • 


Here  X  is  the  free  space  wavelength,  NQ  is  the  maximum  electron 
density  of  the  layer  and  L  is  the  semithickness.  The  situation  is 
illustrated  in  Figure  3.4.  For  unit  input  amplitude  of  the  wave  we 
obtain: 


Figure  3.4  Electron  density  and  dielectric  constant 


plotted  against  z 
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g1  <S,0)  =  2i.(l-S2)  7  ^  •  (^|)T  exp{-ij  -^-((l-S2)1  ♦  1)) 


•  f  AN  (S  i  -6  )  •  Ai|  (-£-)3(5-L-i±T))Ai{(-£-)3(6-L-i!!-T  Cl-S2 ) )  >  d« 


Lc 


Lc 


(3.57) 


where  u>0  is  the  plasma  frequency  at  the  maximum  of  the  layer  where 
N(z)  =  Nc  and  where  Ai(t)  is  the  Airy  integral  defined  by  (14): 


Ai(t) 


}dq 


with  the  path  of  integration  as  shown  in  Figure  3.5. 


(3.58) 


Figure  3.5  Definition  of  path  of  integration 


The  result  (3.57)  agrees  with  that  of  Denisov  (12).  We  first  take  the 
inverse  transform: 
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g,(x,0)  S  (f)3^)1  f"de  Ai{  -  L  -*7)  }  • 

u>0  *  o  Lc  w0 

2  2  u  2  i 

•  f  dS  explik0(Sx-,jL-~(  (1-S2  )^+l)J  }  •  Ai { (— )T(C-L-S^-(1-S2 ) ) }  • 
1  fc>o  Lc  “o2 


•  AN(S,-{)  (3.59) 

The  calculation  of  the  complete  joint  correlation  at  the  reference 
level  is  straight  forward  but  somewhat  laborious.  We  therefore  consi¬ 
der  only  the  case  of  horizontal  stratification.  With  this  assumption 
neglecting  phase  factors,  we  obtain: 


^g j  * ( x , w , 0 )  g(  (x.w+aw.O)^ 


<^(*)3(^)fj  dCd£’  (4NC-{)AN(-5'))a 


•  Ai2{  (— ^— )~(C-L 
Lc2 


2  l 

Ai2l(^r)T(5' 


Lc 


2  1 

“o 


(3.60) 


If  the  waves  penetrate  rather  deeply  into  the  layer,  i.e.  a  distance 

u2 

large  compared  with  the  vertical  scale  m  ,  or  m  <<  L  — -  .  If  this 

“o 

is  fullfilled  we  can  introduce  the  vertical  spatial  correlation  function 
through: 


<^N(-«)AN(-C')^av  =  ^4N(-02^v  exp{-  -y  U-«’)2} 


(3.61) 


where  the  scale  m  is  independent  of  height.  We  obtain: 


38 


^g/(x,u,0)  gj  (x,w+Au,0)^a 


2  2  1 
(Ii£i)7(“)3(]iii1!i)3  f  d?/AN<-02)  Ai2((^_)T(C-L  sL)) 

Z7  0  0  Jo  \  /av  Lc2  Z7 


2  1 

dn  Ai2 ( (—^r) 3 (6  +  n  -  L  ■ u  A  ;^— ) )  exp{  -  n2/m2  } 
Lc2  u,Q2 


(3.62) 


This  double  integral  was  evaluated  by  means  of  an  electronic  digital 
computer  for  a  specific  example.  We  assumed  the  linear  electron  den¬ 
sity  variation  to  extend  over  100  kilometers  and  the  plasma  frequency 
at  the  peak  to  be  5  Mc/s.  Furthermore,  the  vertical  scale  m  was 
allowed  to  assume  the  values  0.1  ,  0.5  ,  1.0  and  10  kilometers.  To 
our  surprise  the  correlation  curves  are  very  insensitive  to  these 
changes  in  vertical  scale  m.  This  is  thought  to  be  caused  by  the  fact 
that  for  irregularities  at  the  reflection  level  these  scales  are  pro¬ 
bably  all  less  than  the  local  wavelength,  and  no  pronounced  effect  of 
scale  change  can  be  expected.  The  total  amount  of  scattered  energy,  on 


Figure  3.6  Correlation  function  for  horizontally  stratified  irregula¬ 


rities  embedded  in  a  linear  electron  density  profile 
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the  other  hand,  was  very  sensitive  to  scale  changes.  The  total  scat¬ 
tered  power  turns  out  to  be  very  nearly  proportional  to  the  vertical 
scale  m.  The  correlation  curve  for  the  above  parameters  are  shown 
in  Figure  3.B.  It  is  seen  to  drop  very  rapidly  to  about  0.6  over 
the  first  10  kc/s  and  from  there  on  the  rate  of  decay  is  very  slow. 
On  the  basis  of  this  model  there  does  not  seem  to  be  much  hope  for 
frequency  correlation  as  a  method  of  studying  the  vertical  scale  of 
the  irregularities.  It  must  be  remembered,  however,  that  the  weak 
scattering  approximation  on  which  the  analysis  is  based  does  break 
down  near  the  reflection  level.  For  frequency  correlation  this  break¬ 
down  might  be  of  comparatively  greater  importance  than  for  space  cor¬ 
relation. 


4  CONCLUSIONS 


In  this  report  we  first  studied  in  rather  general  terms  the  properties 
of  waves  reflected  at  two  fairly  closely  spaced  frequencies  by  a  ran¬ 
dom  medium.  We  described  the  relation  between  the  two  reflected  waves 
by  means  of  a  correlation  function  of  frequency  separation  Au.  In  the 
case  of  an  infinitely  extended  reflector  we  found  that  beyond  a  certain 
critical  distance  from  the  reflector  frequency  correlation  studies  can 
only  provide  the  same  type  of  information  as  spatial  correlation  stu¬ 
dies.  At  smaller  distances  frequency  correlation  can  supply  informa¬ 
tion  in  addition.  For  a  reflecting  medium  of  limited  extent  in  space 
we  were  able  to  show  that  frequency  correlation  tends  to  a  limiting 
value  independent  of  distance  beyond  a  certain  critical  distance  where¬ 
as  the  spatial  correlation  range  would  tend  to  increase  with  distance. 

When  examining  several  ionospheric  models  we  encountered  the  same  dif¬ 
ficulties  as  are  familiar  in  connection  with  the  interpretation  of 
spatial  correlation.  For  deep  phase  modulation  vertical  scale  cannot 
be  deduced  unambiguously  in  the  same  way  as  with  horizontal  scale  in 
diversity  experiments. 
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The  most  important  result  of  the  analysis  presented  above  is  the  sug¬ 
gestion  that  spatial  and  frequency  correlations  should  be  measured 
simultaneously  in  order  to  prevent  angular  spectrum  effects  from  being 
interpreted  as  arising  from  irregular  electron  density  distribution 
along  the  vertical  direction. 
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Part  II 


OBSERVATIONAL  TECHNIQUES  AND  RESULTS 


SUMMARY 


The  present  paper  discusses  four  different  methods  of  study¬ 
ing  tne  correlation  between  signals  reflected  from  the  iono¬ 
sphere  at  closely  spaced  frequencies.  These  methods  are 
applied  to  echoes  from  the  F-layer  at  4  Mc/s.  The  correla¬ 
tion  width  of  the  random  component  is  found  to  be  highly 
variable  between  10  and  100  kc/s  depending  on  the  condi¬ 
tions.  Such  a  frequency  correlation  range  might  be  explain¬ 
ed  by  the  finite  opening  angle  of  the  angular  power  spectrum 
of  the  reflected  waves.  A  significant  specular  component  is 
observed  in  most  of  the  records.  Apparent  vertical  drift 
velocities  as  high  as  40  m/s  are  occationally  recorded, 
but  more  typically  velocities  of  20  -  25  m/s  are  found. 


1  INTRODUCTION 

In  the  present  part  of  this  report  we  discuss  the  principles  of  obser¬ 
vation  of  the  correlation  between  signals  reflected  from  the  iono¬ 
sphere  at  vertical  incidence  at  slightly  different  frequencies.  The 
technical  details  of  the  equipment  used  is  described  in  part  III  of 
this  report.  Here  we  only  deal  with  equipment  questions  insofar  as 
they  are  of  importance  for  the  interpretation  of  the  experimental  data. 

The  results  presented  below  are  only  concerned  with  frequency  correla¬ 
tion.  Due  to  the  somewhat  unexpected  signal  behaviour  described  below 
it  did  not  prove  possible  to  carry  out  simultaneous  spatial  and  fre¬ 
quency  correlation  studies  within  the  time  available.  This  is  parti¬ 
cularly  unfortunate  in  the  light  of  the  theoretical  discussion  of 
part  I  and  in  view  of  the  results  obtained  below. 

In  section  2  we  state  the  observational  problem  arising  in  the  deter¬ 
mination  of  the  correlation  between  the  two  signals.  Several  different 
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techniques  to  study  this  relation  are  then  discussed  in  some  detail.  In 
section  3  we  present  the  observational  results  obtained  by  means  of 
these  techniques,  and  analyse  the  results  with  a  view  to  obtaining  in¬ 
formation  on  some  of  the  properties  of  the  reflecting  medium.  It  should 
be  noted  that  the  results  thus  obtained  are  so  limited  that  the  "typi¬ 
cal"  behaviour  of  the  ionospheric  layers  with  respect  to  frequency  cor¬ 
relation  cannot  be  established.  The  reason  for  this  is  our  inability 
to  pick  observational  periods  at  random  because  of  interference  prob¬ 
lems,  absorption  effects  or  violent  non-stationary  fading.  In  section 
4  we  discuss  the  results  obtained  by  us  and  by  others  in  the  light  of 
the  theories  developed  in  part  I  of  this  report. 


2  METHODS  OF  OBSERVATION 


In  the  theoretical  discussion  of  part  I  of  this  report  we  considered 
the  signals  at  the  two  closely  spaced  frequencies  at  the  receiver  as 
two  random  phasors  rotating  at  different  angular  velocities.  This 
velocity  difference  is  determined  by  the  difference  in  frequency  of  the 
two  sidebands  transmitted.  Suppose  that  this  frequency  difference  is 
removed  after  reception.  We  are  then  left  with  two  phasors  which  are, 
on  the  average  not  rotating  in  any  systematic  manner  with  respect  to 


Figure  2.1  Signal  phasors  at  the  two  frequencies 
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each  other  as  long  as  the  reflecting  medium  remains  statistically 
stationary.  The  two  phasors  are  represented  either  by  their  x,y 
components  or  by  their  polar  coordinates  A,*  as  shown  in 
Figure  2.1'. 

The  statistical  relationship  between  the  two  phasors  can  be  described 
by  means  of  the  quantities: 

(xix2)av  •  (yiy2)av 

(2.1) 

<x.y2)av  •  (x*y.)av 

Because  of  symmetry  we  shall  assume  that: 

(*i*2)av  *  (yi^av 

and  (2.2) 

(xiy2)av  =  -(x2yi)av 

In  terms  of  the  polar  coordinate  representation  we  have: 

(xlxl)av  2  j(a,A2(cos(»j-*,)  +  cos(*l+»2)))av  = 

2  i<A,A2  cos(*2-»1))av 

(yiy2)av  2  J(a,A2(cos(#2-*1)  -  cos(*,+*2)))av  = 

2  KA‘Aj  c°s(*i-*l))av 

\xiy2)av  =  i(AiA2(sin(*2-*,)  +  sin(*, +  *2 ) ))av  = 

2  KAiA2  sin(*2-*:>)av 

(x2yi)av  s  i(+A,A2(sin(*,+*2)  +  sin(*,-#2) ))ay  * 

=  -i(A1A2sin(#2-*1))av 


(2.3) 
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This  is  related  to  the  complex  representation  of  part  I  through: 


g,  =  —  A,  exp{  i  } 


5  =  A  expl  i  4  } 

2  /J  2  *2  ' 


(2.4) 


Using  the  definition  that  ^gtg2^>av  =  Rg  We  find  that 

(xix2)av  1  (y^av  =  Re<Rg> 


and 


(2.5) 


(V2)av  =  -<X2y:Xv  =  Im(Rg 


In  the  autocorrelation  analysis  of  a  noise  signal  the  two  phasors  which 
are  correlated  are  taken  from  a  single  signal  and  the  average  phase 
difference  of  the  two  phasors  is  zero.  In  the  case  considered  here  the 
two  phasors  are  derived  from  different  signals  which  have  been  reflec¬ 
ted  from  the  ionosphere  at  different  frequencies.  This  means  that  a 
systematic  phase  difference  does  exist  between  the  two  and  that  this 
difference  depends  strongly  on  the  modulation  frequency  if.  This  can 
be  seen  most  easily  by  means  of  the  following  model:  Suppose  the  iono¬ 
sphere  is  a  plane,  perfectly  reflecting  surface  at  height  H.  In  this 
case  the  phase  difference  becomes: 

4* 

♦2  -  ♦,  =  4*  *  —  H  .  if  (2.6) 

Therefore,  even  with  relatively  small  large-scale  variations  in  height 
the  two  phasors  will  rotate  with  respect  to  each  other.  This  in  turn 
means  that  the  two  phasors  to  be  correlated  are  not  related  to  each 
other  in  the  same  simple  way  the  two  phasors  occurring  in  the  auto¬ 
correlation  analysis  of  gaussian  noise.  The  rate  of  rotation  A*  is 
related  to  the  apparent  vertical  velocity  through: 
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(2.7) 


Because  of  this  phase  uncertainty  the  distinction  between  real  and 
imaginary  parts  of  R^  looses  its  significance,  and  it  is  the  magni¬ 
tude  j  Rg |  which  really  matters. 

In  the  discussion  so  far  we  have  imagined  the  random  signal  to  have  pro¬ 
perties  like  those  of  band-pass-limited  gaussion  noise.  As  will  be 
shown  below  there  is  in  many  cases  a  systematic  "specular"  component 
present  in  addition  to  the  random  noiselike  component  described  above. 
Such  a  specular  component  could  arise  through  reflection  at  a  shallow 
phase  modulating  reflector  as  we  saw  in  part  I  of  the  present  report. 

The  presence  of  such  a  specular  component  has  important  implications 
on  the  determination  of  |Rg|  for  the  random  component  from  post-de¬ 
tection  data.  The  quantity  which  can  be  determined  directly  by  the 
methods  described  below  is  the  post-detection  correlation  coefficient 
defined  by: 


ra  =  (GlG2)av  ~  ^,G‘.)av(gj}av  (2.8) 

-(Gi>av>%  -<G2)av>2)av 

and  the  relationship  between  the  pre-  and  the  postdetection  correlation 
will  depend  to  some  extent  on  the  ratio  of  the  power  in  the  "specular" 
and  the  random  components: 


where  Pj  and  P2  are  the  amplitudes  of  the  specular  components  at 
the  two  frequencies,  and  Aj  and  A2  those  of  the  random  part.  G, 
and  G2  are  the  resultant  amplitudes. 
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Yet  another  difficulty  arises  in  the  correlation  of  signals  at  spaced 
frequencies  when  a  "specular"  component  is  present.  The  specular  com¬ 
ponents  might  have  arbitrary  phase  relationships  depending  on  the  dis¬ 
tance  to  the  reflector  as  just  explained.  If  the  random  components  do 
not  arise  through  scattering  or  reflection  at  the  same  distance  as  the 
specular  components  the  relative  phases  of  the  random  and  the  specular 
components  in  the  two  sidebands  might  be  different  and  serious  diffi¬ 
culties  might  arise  in  the  determination  of  the  relation  between  pre- 
and  post-detection  correlations.  Consider  as  an  example  the  particu¬ 
lar  case  of  the  two  random  components  being  completely  correlated 
(identical)  and  much  smaller  than  the  "specular"  components.  When 
the  two  "specular"  signals  are  in  phase,  the  post-detection  correlation 
will  be  unity.  When  they  are  90  or  270°  out  of  phase  :  0 
and  when  the  phase  difference  is  180°  R  *  -1  .  Clearly,  there- 


Figure  2.2  Relation  between  pre-  and  postdetection  correlation  for 
zero  phase  shift  with  a  =  specular/random  as  parameter 
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fore,  one  must  be  very  careful  in  drawing  conclusions  about 


The  relation  between 


and  has  been  computed  by 


A . . .  '  g'  A  -  r 

K  Bl^takjar  <1)  for  the  case  of  zero  phase  difference  between  spe 

cular  components  with  o  as  parameter.  The  curves  are  shown  in 


Figure  2.2. 


Because  of  the  importance  of  knowing  whether  a  specular  component  is 
present  or  not  the  results  were  checked  for  a  specular  component  in 
the  following  way:  In  order  to  compute  from  our  data  it  was  ne¬ 
cessary  to  find  (GlG2)av  »  (Gl)av  >  (Gf)av  »  {Gl)av  and  .  From 

simple  application  of  the  theory  of  gaussian  noise  with  a  "specular" 
component  (2)  we  can  find  a  relation  between  the  quantity 
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-AWav -  (2>10) 

(G?)av(G2)av 


and  Ra  provided  (s,)av  =  (s2)av  and  (Gl)av  =  (G22)av 
tion  depends  on  a  as  defined  in  (2.9): 


This  rela- 


R0  ■  KA*(1-RA,rsfr7|(1‘“)Io(2)*lIi(!))!  <2-n> 


where  IQ  and  I(  are  Bessel  functions  with  imaginary  arguments. 
This  relation  is  plotted  in  Figure  2.3  with  a  as  parameter. 


2.1  Direct  amplitude  correlation 

The  most  obvious  method  to  determine  the  correlation  between  signals 
reflected  at  slightly  different  frequencies  is  to  correlate  their  amp¬ 
litudes  directly.  As  carried  out  by  us  the  two  sidebands  were  detected 


Figure  2.4  Block  diagram  of  direct  amplitude  correlation 
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in  separate,  independent  receivers  tuned  one  to  each  of  the  sidebands. 
The  detected  outputs  were  fed  to  a  dual  channel  pen  recorder.  Corre¬ 
lation  coefficients  were  derived  by  reading  off  the  amplitude  values 
and  correlating  the  two  sets  of  readings  by  means  of  a  desk  calculator. 
The  system  is  shown  in  Figure  2.4.  The  system  works  satisfactorily 
but  suffers  from  the  drawback  that  no  phase  information  is  available. 
When  used  in  connection  with  calculations  by  hand  the  method  is  quite 
laborious.  Direct  amplitude  correlation  was  only  used  by  us  in  order 
to  check  other  methods. 


2 . 2  Correlation  by  sum  and  difference  of  signals 

Because  of  the  desirability  to  keep  track  of  the  phase  relationship  of 
the  two  sidebands  we  first  tried  a  method  of  observation  based  on  prin¬ 
ciples  related  to  interferometry.  On  reception  of  the  two  sidebands 
we  first  removed  the  systematic  frequency  difference  between  the  two 
signals  coherently  as  shown  in  principle  in  Figure  2.5.  The  difference 
frequency  at  the  transmitter,  af  ,  is  derived  from  the  difference  of 
the  frequencies  of  the  two  receiver  local  oscillators  L01  and  L02. 
The  two  IF  signals  were  then  applied  to  a  circuit  shown  in 
Figure  2.6.  This  circuit  forms  the  sum  and  the  difference  of  the  two 
signals  with  zero  phase  shift  and  with  90°  phase  shift  applied  to 
one  of  the  signals  before  adding  and  subtracting.  Further  details  of 
the  equipment  can  be  found  in  part  III  of  the  report. 

By  combining  the  signal  phasors  gj  and  g2  as  shown  in  Figure  2.6 
we  obtain  after  detection  in  envelope  detectors: 


(6? 


+  g: 


2GjGj  cos(e2 


-  e,»i 


(2.12a) 


and 


50 


Figure  2.5  Block  diagram  showing  the  principle  of  bringing  the 
two  sidebands  coherently  to  the  same  frequency 


(G*  +  G j  ‘  2Gj G2  sin(e2  -  e,))i 


(2.12b) 


(We  use  G  and  e  here  to  indicate  total  signal,  possibly  consisting 
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of  both  random  and  specular  components).  By  adjusting  the  gain  of  the 
receivers  or  the  relative  strength  of  the  sidebands  transmitted  one 
can  make  the  average  values  of  Gj  and  G2  equal.  This  means  that 
we  can  find,  by  proper  averaging  the  following  quantities. 


If  the  signal  were  of  a  gaussian  type  with  no  specular  component  these 
expressions  would  give  the  real  and  imaginary  parts  of  the  complex  cor¬ 
relation  coefficient  Rg.  With  a  "specular"  component  we  would  get 
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Figure  2.7  Example  of  a  record  of  S  and  of  the  amplitude 
G  of  one  of  the  signals 


into  trouble  unless  the  relative  phase  of  the  two  "specular"  signals 
were  known. 

In  practice  it  is  not  so  simple.  Due  to  large  scale  changes  in  height 
the  phase  of  the  two  phasors  will  change  with  respect  to  each  other  and 
we  always  find  that  ^S2^av  1  Saving  zero  correlation.  Figure 

2.7  shows  an  example  of  a  record  of  S  and  of  the  amplitude  G  of  one 
of  the  signals.  We  see  that  the  record  of  G  shows  very  little  fading. 
S  on  the  other  hand  during  the  period  shown  in  Tigure  2.7  goes  through 
more  than  one  complete  fading  cycle,  which  can  only  be  explained  in 
terms  of  rotation  of  one  of  the  phasors  with  respect  to  the  other  one. 
For  this  particular  record  as  many  as  40  complete  revolutions  were 
observed  during  the  course  of  one  hour.  This  record  was  taken  on 
October  5  1962  at  10-11  hrs.  LT  at  a  frequency  of  4  Mc/s  and  with 
4f  =  31  kc/s.  The  phase  change  versus  time  is  shown  in  Figure  2.8,  on 
the  assumption  of  rotation  in  the  same  direction  all  the  time.  There 
is  an  ambiguity  in  the  direction  of  rotation  in  this  type  of  record. 
Round  the  points  A  and  B  in  Figure  2.8  it  might  be  that  the  direc¬ 
tion  of  rotation  reverses.  -  On  the  assumption  of  rotation  in  the 
same  direction  all  the  time  the  vertical  velocity  comes  out  to  be  of 
the  order  of  25-30  m/s  which  is  rather  high,  but  not  impossibly  high 
when  compared  with  observed  vertical  velocities  (3,4).  This  type  of 
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behaviour  was  more  or  less  pronounced  in  all  the  records  made  by  this 
method.  Because  of  this  phase  rotation  the  evaluation  of  the  magni¬ 
tude  of  the  correlation  coefficient  would  require  the  simultaneous  re¬ 
cording  of  S,D,S  and  D.  Due  to  difficulties  in  keeping  the  gain  of 
four  channels  exactly  equal  and  due  to  a  rather  tedious  analysing 
technique,  this  method  was  abandoned  in  the  present  series  of  experi¬ 
ments  . 


Figure  2.8  Revolution  of  ohasors  as  a  function  of  time 


2 . 3  Frequency  sweep  technique 

In  the  expressions  (2.12a)  or  (2.12b)  the  phase  difference 

49  =  6j  -  6j  is  determined  by  the  distance  to  the  reflecting  medium 

roughly  according  to  formula  (2.6): 

49  *  —  •  H  •  4f 

c 
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For  S  or  D  we  therefore  obtain: 


S 


D 


(<3?  * 


G2  *  2GjG2  COS 


(T«  • 


(2.14) 


Suppose  now  that  we  let  if  vary  with  time,  in  a  linear  fashion,  say. 
The  result  will  be  that  a  modulation  will  be  imposed  on  S  and  that 
the  depth  of  this  modulation  will  depend  on  the  product  GjG2.  This 
incidentally  under  quiet  conditions  can  be  used  to  determine  quite  ac¬ 
curately  the  apparent  height  of  reflection  of  the  signal  (h1). 

Figure  2.9  shows  an  example  of  part  of  such  a  frequency  sweep.  The 
upper  channel  indicates  the  frequency  difference  Af  and  the  lower 
channel  the  sum  S.  The  irregular  nature  of  the  record  is  caused  by 
the  manual  frequency  variation.  By  counting  26  complete  relative  re¬ 
volutions  the  virtual  height  was  determined  to  be  228  km  at  the  fre¬ 
quency  used,  fQ  =  4  Mc/s. 


Figure  2.9  Determination  of  virtual  height  bv  frequency 
change  technique 


By  changing  the  modulation  frequency  Af  periodically  by  about 
1500  c/s  round  some  modulation  frequency  and  by  noting  the  values  of 
the  maxima  and  the  minima  it  is  possible  to  determine  G2  ♦  G2  and 
|G,  -  Gj|  .  By  forming  the  quantity 
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(G,  ♦  G2 ) 2  -  | G[  -  G2 | 2 

-  *  R0  (2.15) 

(Gj  ♦  G2 ) 2  ♦  | Gj  -  Gj | 2 

(see  (2.10))  the  post-detection  correlation  function  can  be  found 
from  (2.8)  by  also  evaluating 

<Si)av  *  «>av  »  KG‘+G^)av  <2-16) 

This  method  gives  no  more  information  than  the  direct  correlation  of 
individually  received  signals,  but  it  turns  out  to  be  more  convenient 
for  manual  data  reduction  because  the  upper  and  lower  envelopes  of  the 
signal  modulation  can  be  treated  individually.  An  example  of  such  a 
record  is  shown  in  Figure  2.10.  About  half  of  the  analysed  frequency 
correlation  data  were  obtained  by  this  technique,  see  table, 

Appendix  I . 


Figure  2.10  Sample  of  frequency  sweep  record 


2.4  Phase  changing  technique 

The  sum  and  difference  technique  described  in  section  2.2  can  be  modi¬ 
fied  somewhat  as  shown  in  Figure  2.11,  The  only  difference  is  the  ad¬ 
dition  of  a  continuous  phase  changer  in  one  of  the  signal  paths.  By 
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rotating  this  phase  changer  at  a  constant  rate  one  of  the  signal  pha- 
sors  can  be  made  to  rotate  with  respect  to  the  other.  By  recording  the 
sum  S  as  before,  it  is  possible  to  observe  the  modulation  caused  by 
this  rotation,  and  by  reading  off  the  maxima  and  the  minima  as  in  the 
frequency  changing  method  the  amplitude  correlation  function  can  be 
found.  An  example  of  such  a  record  is  shown  in  Figure  2.12.  On  the 
other  channel  there  are  spikes  which  serve  as  references  to  indicate 
the  position  of  the  phase  changer.  There  is  one  such  spike  per  ^ 
radians  of  phase  change,  see  Figure  3.M.  Phase  information  can  be  de¬ 
rived  by  considering  the  relative  position  of  the  modulation  pattern 
of  the  signal  pattern  and  the  reference  spikes.  For  the  analysis  of 
such  records  we  had  available  an  analogue  to  digital  converter  and  a 
tape-recorder  to  store  the  digital  data.  The  digital  data  were  then 
analysed  in  a  computer.  The  following  quantities  were  derived: 


(2.17) 
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and 

(2.18) 

Furthermore,  the  relative  phase  46  of  the  reference  spikes  and  the 
modulation  every  ten  seconds,  and  finally  an  expression  for  the  random 
phase  fluctuations  viz. 


This  latter  quantity  is  rather  difficult  to  interpret  in  that  it  is 
jointly  determined  by  time  variation  (time  between  sampling  points 
about  1  sec)  and  variation  of  correlation  with  frequency. 


Observations  of  the  correlation  between  the  signal  reflected  at  spaced 
frequencies  at  vertical  incidence  from  the  ionosphere  were  begun  during 
September  1962  and  continued  into  the  first  weeks  of  January  1963 
at  Kjeller,  near  Oslo.  All  the  methods  described  above  except  the 
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straight  sum  and  difference  procedure  were  used  for  qualitative  eva¬ 
luation  of  the  correlation.  In  connection  with  the  phase  changing 
method  electronic  sampling  and  digital  computer  analysis  was  used. 

Other  records  were  analysed  by  hand. 

The  frequency  used  was  4  Mc/s  throughout  the  experiment,  this  being 
the  frequency  available  which  was  least  susceptible  to  interference. 

The  other  two  frequencies  were  2.7  and  5.4  Mc/s.  During  the  obser¬ 
vational  period  the  critical  frequency  of  the  F-layer  was  in  the  re¬ 
gion  5  -  6.5  Mc/s  and  E-layer  frequencies  less  than  4  Mc/s.  Our 
results  -  taken  during  the  daytime  -  therefore  pertain  to  the  F-layer. 
Several  attempts  were  made  at  observing  during  the  nighttime,  but  this 
proved  impossible  because  of  severe  interference  problems.  All  the  re¬ 
sults  presented  below  therefore  apply  for  F  layer  day-time  conditions. 

For  all  the  successful  records  the  quantities  Rq  and  RA  were  calcu¬ 
lated  as  presented  in  the  Table  in  Appendix  I  at  the  end  of  this  part 
of  the  report.  From  RQ  and  RA  diagrams  showing  RQ  as  a  function 
of  Ra  were  plotted  for  all  the  records  and  the  strength  of  the  spe¬ 
cular  component  was  estimated  from  a  comparison  with  Figure  2.3.  This 
is  discussed  in  section  3.1  below.  Next  those  data  which  could  be  con¬ 
verted  into  predetection  correlation  were  picked  out  for  special  study 
as  described  in  section  3.2.  Finally,  in  section  3.3,  some  results 
from  the  study  of  phase  variation  of  the  signal  are  briefly  considered. 


3.1  Deductions  about  specular  component 

Figure  3.1  shows  two  examples  of  the  evaluation  of  the  ratio  of  specu¬ 
lar  power  to  random  power  in  the  returned  signal.  The  lines  drawn  in 
full  correspond  to  the  best  fit  with  the  observational  data.  Some  of 
the  curves  of  Figure  2.3  are  shown  dotted  in  the  same  diagram.  In  most 
of  the  records  a  distinct  specular  component  was  present  and  in  view 
of  the  possibility  of  the  specular  components  to  be  out  of  phase  at  the 
two  sideband  frequencies  we  did  not  dare  to  make  any  deductions  about 
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Figure  3.1  Evaluation  of  specular  to  random  power  ratio 


for  records  10  and  16 
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the  predetection  correlations  for  the  random  component  except  for  very 
small  values  of  a.  From  the  table  in  Appendix  I  it  is  apparent  that 
the  records  with  phase  changing  (No  13-20)  nearly  all  have  a  strong 
specular  component.  This  is  thought  to  be  due  to  our  tendency  to  make 
these  observations  under  conditions  of  fairly  slow  fading,  the  reason 
for  this  being  that  the  imposed  rate  of  phase  rotation  should  prefer¬ 
ably  be  much  higher  than  the  rate  of  fading  of  the  signal.  This  was 
unfortunately  only  realised  after  the  recordings  were  completed. 


3.2  Correlation  coefficients 


From  the  table  of  records.  Appendix  I,  it  is  apparent  that  only  the 
following  records  can  be  used  for  the  evaluation  of  the  correlation 
coefficient  of  the  random  reflected  component  at  spaced  frequencies: 
records  number  4,  6,  9,  10,  12  and  19.  The  consistent  correlation 
coefficients  from  these  records  were  plotted  in  Figure  3.2.  No  attempt 
was  made  to  combine  the  data  into  one  curve  showing  correlation  against 
frequency  separation  because  it  was  realized  that  condition  under  which 
the  data  were  taken  could  well  be  as  different  as  is  apparent  from  the 
observational  results.  For  each  of  the  useful  records  an  attempt  was 
made  to  draw  a  curve  of  correlation  versus  frequency  separation.  It 
was  considered  sufficiently  accurate  to  draw  this  curve  according  to 
the  judgement  of  the  authors  rather  than  according  to  some  least  mean 
square  error  criterion  or  the  like.  There  are  particularly  two  points 
at  Af  -  66  kc/s  in  record  9  which  do  not  fit  the  curves  at  all. 
Otherwise  the  points  fit  the  curves  quite  well.  In  order  to  see  whether 
the  spread  in  the  data  could  be  explained  from  statistical  errors  caused 
by  the  finite  number  of  sampling  points  used  in  the  evaluation  of  the 
correlation  coefficient  we  calculated  the  upper  and  lower  quartiles  for 
the  correlation  coefficient  on  the  assumptions  that  a)  the  two  sig¬ 
nals  correlated  originate  in  a  bivariate  normal  population  with  corre¬ 
lation  coefficient  p  ,  b)  the  correlation  coefficient  p  is  identi¬ 
cal  with  the  one  shown  by  the  dotted  curves  in  the  diagrams  and  c) 
that  the  number  of  independent  pairs  of  samples  giving  rise  to  each 
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point  in  the  diagram  is  20,  which  roughly  corresponds  to  the  typical 
number  of  fading  cycles  in  a  record.  It  is  seen  that  most  of  the 
points  fall  inside  this  interquartile  range.  The  random  spread  in  the 
data  of  the  different  diagrams  can  therefore  well  be  explained  as  pure¬ 
ly  statistical  fluctuations  caused  by  the  finite  number  of  sampling 
points.  Details  of  the  evaluation  of  the  interquartile  range  can  be 
found  in  Appendix  II.  The  only  other  source  of  error  which  is  consi¬ 
dered  important  in  the  evaluation  of  correlation  coefficients  is  the 
possibility  of  unequal  gain  in  the  two  channels.  As  shown  in 
Appendix  III  this  could  lead  to  an  error  in  ,  a Rft  given  by: 


AR* 


Ra  -  R1 

A  A 


4-ff 


{»  ♦  R. C  8 — w )  } 


1  +  i  + 


8-ii 

4-w 


(3.1) 


Rg 


Figure  3.3  Systematic  error  in  R^  due  to  unequal  channel  gain 
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on  the  assumption  of  no  specular  component.  Here  R^  is  the  apparent 
and  the  "real"  envelope  correlation  and  i  is  the  excess  gain  of 
one  of  the  channels.  The  corresponding  error  in  the  predetection  cor¬ 
relation  is  approximately: 


AR 


g 


(3.2) 


This  function  is  shown  in  Figure  3.3  for  two  different  values  of  6. 

In  our  experiment  6  is  generally  smaller  than  0.1.  The  systematic 
error  due  to  unequal  gain  is  therefore  much  less  than  the  random 
errors . 

The  errors  introduced  through  noise  and  interference  are  thought  to  be 
unimportant  because  records  were  only  taken  when  the  signal  to  noise 
ratio  was  good,  (better  than  10/1  in  amplitude).  It  is  therefore 
concluded  that  the  results  in  Figure  3.2  represent  the  real  correlation 
coefficient  subject  only  to  the  random  uncertainties  indicated. 


3 . 3  Phase  Variation 


As  explained  in  section  2.4  the  phase  changing  method  should  give  an 
opportunity  for  tracking  the  phase  difference,  and  the  phase  was  ac¬ 
tually  printed  out  for  every  ten  seconds.  Upon  closer  examination  of 
the  record  it  appeared  that  the  phase  changes  printed  out  from  the 
sampler-computer  combination  were  rather  too  large  -  often  leading  to 
apparent  vertical  velocities  of  the  order  of  150  -  200  m/sec.  A 
thorough  examination  of  the  record  samples  revealed  that  once  in  a 
while  a  sample  was  missing  or  an  extra  sample  was  inserted  due  to  some 
sort  of  spurious  noise  in  the  sampling  device.  This  has  lead  to  phase 
jumps  of  90°  here  and  there  in  the  analysed  results.  For  this  rea¬ 
son  we  have  rejected  all  the  phase  information  in  the  results  analysed 
automatically.  As  an  example  of  the  type  of  information  which  is 
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Figure  3.4  Sample  from  phase  changing  record  showinR  rapid 
phase  change 


Figure  3.5  Relative  phase  change  versus  time  deduced  from 


the  record  shown  in  Figure  3.4 


Apparent  height  change  (km) 
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inherent  in  the  phase  changing  data  we  have  analysed  a  few  records  by 
hand  in  order  to  evaluate  possible  systematic  apparent  vertical  velo¬ 
cities.  An  example  of  a  particularly  rapid  phase  variation  is  shown  in 
Figure  3.4.  Part  of  the  corresponding  record  can  be  seen  in  Figure  3.5. 
The  average  velocity  over  the  period  of  this  record  is  48  m/sec  direc¬ 
ted  upwards.  This  is  exceptionally  high.  More  often  we  find  veloci¬ 
ties  of  the  order  of  25  m/sec.  It  is  thought  that  part  of  the  high 
apparent  velocities  might  be  related  in  some  way  to  the  presence  of  only 
a  few  components  in  the  angular  spectrum  and  to  an  angular  flicker  of 
the  reflected  signal  which,  in  itself  quite  well  could  give  rise  to  high 
apparent  vertical  velocities  over  short  time  intervals.  Obviously  the 
simultaneous  observation  of  frequency  correlation  and  direction  of  ar¬ 
rival  would  give  a  very  detailed  picture  of  the  reflection  mechanism. 


4  DISCUSSION  AND  CONCLUSIONS 


In  the  present  part  of  the  report  we  were  first  able  to  show  that  in 
most  of  the  cases  considered  here,  viz  day-time  reflections  from  the 
F-layer,  there  is  an  appreciable  specular  component  present.  The  pre¬ 
sence  of  this  component  makes  the  interpretation  of  the  observations 
difficult  because  there  might  be  a  different  phase  relationship  be¬ 
tween  the  specular  components  than  between  the  random  component  we  are 
interested  in. 

When  no  appreciable  specular  component  was  present  and  hence  when  re¬ 
liable  interpretation  of  the  data  could  be  carried  out  the  correlation 
appeared  to  drop  to  about  0.5  over  a  frequency  separation  range  some¬ 
where  between  about  10  kc/s  and  100  kc/s.  This  result  appears  to  be 
in  rather  good  agreement  with  the  results  of  Gusev  and  Mirkotan  (3). 
Before  conclusions  are  drawn  about  vertical  scales  the  above  values 
should  be  interpreted  in  terms  of  angular  spectra  of  waves  as  explained 
in  part  I  of  this  report.  It  turns  out  that  the  10  kc/s  correlation 
width  can  be  explained  by  a  8m  of  23°  and  the  100  kc/s  by  7.3°. 
If,  in  addition,  it  is  remembered  that  our  theoretical  considerations 


66 


are  based  on  a  one-dimensional  analysis,  we  may  conclude  that  the  ac¬ 
tually  required  values  of  6m  are  even  smaller  than  these  because  of 
the  two-dimensional  nature  of  the  actual  reflectors.  These  values  come 

very  close  to  the  6  which  have  been  observed  by  others  for  waves  re- 
in 

fleeted  from  the  F-layer  under  similar  conditions  (4,5).  It  therefore 
appears  to  us  that  no  conclusions  about  vertical  scales  should  be  drawn 
before  simultaneous  observations  of  angular  spectrum  and  frequency  cor¬ 
relation  can  be  carried  out. 

Finally,  our  observations  of  the  very  large  fluctuations  and  changes 
in  the  virtual  height  of  reflection  of  radio  waves  from  the  ionosphere 
probably  mean  that  the  ionosphere  must  be  regarded  as  a  large  scale 
deeply  corrugated  layer  with  smaller  scale  irregular  structure  super¬ 
imposed. 
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APPENDIX  II 


In  order  to  evaluate  the  upper  and  lower  quartiles  for  the  correlation 

coefficient  we  follow  the  method  described  by  Fischer  (1).  We  assume 

that  the  correlation  coefficient  R  is  a  random  variable  with  mean 

8 

value  p  .  It  is  furthermore  assumed  that  Rg  arises  from  a  finite 
sample  of  pairs  of  values  from  a  bivariate  normal  distribution  with 
correlation  coefficient  p  .  According  to  Fischer  the  quantity 


z 


i 


In 


1  +  R 
1  -  R 


£ 

g 


is  approximately  normally  distributed  about  a  mean  value 


(z)av  *  i in  rHr 

with  a  variance  —~r  where  n  is  the  number  of  pairs  of  values  used 
n-3 

for  the  evaluation  of  R  ,  In  our  case  n  *  20  ,  and  the  variance 
g 

M  *  0.0588.  The  standard  deviation  o  *  0.24.  Now,  upper  and  lower 
quartiles  deviate  from  the  mean  value  by  0.6745  o  *  0,16.  In  the  fol¬ 
lowing  table  we  give  the  upper  and  lower  quartiles  for  a  set  of  diffe¬ 
rent  values  of  the  correlation  coefficient  P  of  the  mother  populations: 


Table  A2 


p 

z  av 

Z 

upper 

quart 

Rg 

upper 

quart 

z 

lower 

quart 

Rg 

lower 

quart 

0.9 

1.47 

1.63 

0.93 

1.31 

0.86 

0.8 

1.10 

1.26 

0.85 

0.94 

0.74 

0.7 

0.87 

1.03 

0.77 

0.71 

0.61 

0.6 

0.69 

0.85 

0.69 

0.53 

0.49 

0.5 

0.55 

0.71 

0.61 

0.39 

0.37 

0.4 

0.42 

0.58 

0.52 

0.26 

0.25 

0.3 

0.31 

0.47 

0.44 

0.15 

0.15 

0.2 

0.20 

0.36 

0.35 

0.04 

0.04 

0.1 

0.10 

0.26 

0.25 

-0.06 

-0.06 

Reference 


(1)  Fischer,  R  A 


"Statistical  Methods  for  Research  Wor¬ 
kers”.  Oliver  and  Boyd,  Edinburgh  (1938) 
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APPENDIX  III 


Assume  that,  due  to  unequal  gain  in  the  system,  the  two  observed  signal 
amplitudes  a(  and  a2  are  given  by 


a2  =  A2<1  ♦  5) 

where 

(Ai)av  1  (A  2 )a  v  =  (A)av 
(A?)av  =  <A0av  * 

and  where  S  is  the  fractional  gain  difference.  If  we  evaluate  the 
envelope  correlation  according  to  (2.18)  the  apparent  envelope  corre¬ 
lation  becomes: 


^(A|+A2(l  +  «))2  -  (Ai-A2Cl  +  <S))2^av  -  l^Ai+A2(l  +  «)^av)2 
((A,+A2<l  +  «>)2  +  (A,-A2(l  +  6))2^av  -  l(A,+A2(l+«)^av)2 


P.A(l  +  «) 


1  ♦  S  * 


■■  K>2).v  - 

*  <*>  - 


Using  the  following  expressions  for 


and 


2 


<Ai)av  *  2  *o 
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KA)av)2  =  f  *o 

and  making  use  of  the  approximate  relation 


the  result  (3.1)  is  obtained. 
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Part  III 


DESCRIPTION  OF  EXPERIMENTAL  EQUIPMENT 


SUMMARY 


A  block  diagram  of  the  complete  experimental  equipment  is 
presented  and  the  principle  of  operation  are  explained. 

The  equipment  is  composed  of  three  main  parts,  and  the  de¬ 
scription  of  these  are  given  with  reference  to  detailed 
circuit  diagrams.  The  procedure  for  adjusting  the  equip¬ 
ment  is  briefly  described. 


1  INTRODUCTION 


This  part  of  the  report  describes  the  experimental  equipment  used  in 
the  measurement  of  correlation  between  radio  signals  reflected  from 
the  ionosphere  at  adjacent  frequencies. 

Most  of  the  equipment  is  of  conventional  design,  and  only  those  fea¬ 
tures  that  are  specific  to  the  experiment  are  dealt  with  in  some  de¬ 
tail.  The  description  will  be  given  with  reference  to  detailed  circuit 
diagrams  and  block  diagrams  that  depict  the  principles  of  operation. 


2  BLOCK  DIAGRAM 


In  this  section  a  brief  description  of  the  block  diagram  of  the  com¬ 
plete  experimental  system  (Figure  2.1)  will  be  given. 

As  indicated  in  the  diagram,  the  equipment  is  subdivided  into  three 
main  parts  i.e.:  the  pulse  transmitter,  the  receiver  equipment,  the 
detector  and  recording  system. 
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The  pulse  transmitter  is  located  in  a  separate  house  to  shield  the  re¬ 
ceivers  from  the  transmitted  pulse.  The  modulation  frequency  f^  is 
sent  via  cable  to  the  transmitter.  The  carrier  frequency  is  suppres¬ 
sed  and  the  transmitted  pulses  therefore  contain  only  two  adjacent  fre¬ 
quencies  fQ  *  f^. 

From  the  receiving  aerial  the  signal  passes  through  a  gate  to  the  re¬ 
ceivers  where  the  two  sidebands  present  in  the  echo  pulses  are  separa¬ 
ted.  To  preserve  the  phase  relationship  between  the  two  echo  signals 
in  the  IF  channel,  the  mixer  frequencies  must  be  coherent.  The  sub¬ 
unit  on  the  upper  left  side  of  the  diagram  generates  these  coherent 
frequencies  and  also  the  modulation  frequency  fL.  A  second  gate  is 
inserted  into  the  IF  channel.  Both  this  gate  and  the  antenna  gate 
are  controlled  by  the  pulse  generator,  the  open  period  coinciding 
with  the  reflected  pulse. 

In  the  detector  and  recording  system  the  signals  in  the  IF  channels 
are  combined,  detected  and  finally  recorded.  To  measure  the  phase 
between  the  two  signals  they  are  combined  with  different  phase  rela¬ 
tionships  by  the  phase  changer  and  the  addition  circuit.  The  detected 
video  pulses  are  passed  to  a  box  car  detector,  which  is  a  special  sort 
of  peak  voltmeter.  This  peak  detected  signal  is  recorded,  in  analog 
form  on  a  Brush  recorder  and  in  digital  form  on  a  tape  recorder  equip¬ 
ped  with  an  analog-digital  converter. 


3  DESCRIPTION  OF  THE  PULSE  TRANSMITTER 


This  section  describes  the  pulse  transmitter  and  its  auxiliary  ampli¬ 
fiers  with  reference  to  three  circuit  diagrams  (Diagrams  3.1  -  3.3). 
The  aerial  system  is  also  described. 
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3.1  Block  diagram 

Figure  3.1  gives  a  more  detailed  block  diagram  of  the  pulse  transmitter. 
The  carrier  is  generated  in  a  crystal  oscillator  and  fed  to  a  balanced 
modulator  through  a  gated  amplifier.  The  modulator  output  is  amplified 
in  two  stages  buffer  and  power  amplifier.  The  latter  is  gated  with  the 
same  pulsed  signal  as  the  hF-amplif ier. 


Figure  3.1  Block  diagram  of  pulse  transmitter 


3 . 2  HF  oscillator 


The  HF  oscillator  (Diagram  3.1)  can  generate  three  fixed  frequencies. 
It  is  a  type  of  electron  coupled  oscillator  with  crystals  as  the  fre¬ 
quency  sensitive  elements.  The  HF  signal  is  amplified  in  one  stage  and 
the  output  is  taken  from  a  cathode  follower. 
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3 . 3  Pulse  generator 

A  monostable  multivibrator  is  triggered  from  the  A  C  mains  (SO  c/s) 
with  adjustable  phase.  (Diagram  3.1).  The  pulse  duration  is  variable 
between  0.5  and  4  msec.  The  output  pulses  are  amplified  and  brought 
out  via  a  cathode  follower. 

3.4  LF  amplifier  and  phase  inverter 

The  input  to  the  LF  amplifier  (Diagram  3.2,  upper  left  side)  is  the 
modulation  frequency  f^.  A  long  tailed  pair  coupling  is  used  for 
amplification  and  phase  splitting. 

3. 5  HF  amplifier 

This  amplifier  (Diagram  3.2)  consists  of  a  long  tailed  pair  which  is 
gated  by  current  gating  to  the  cathodes.  The  output  is  taken  through 
a  transformer  to  prevent  the  gating  pulse  from  reaching  the  modulator. 
The  output  circuit  is  heavily  damped  to  obtain  a  good  step  response. 

3. 6  Balanced  modulator 

The  balanced  modulator  (Diagram  3.2)  consists  of  two  valves.  The 
carrier  is  applied  to  the  control  grids  with  equal  phases  and  the  modu¬ 
lating  signal  to  the  screen  grids  with  opposite  phases.  The  anode  load 
is  centertapped  and  transformer  coupled  to  the  next  stage. 

3. 7  Buffer  and  power  amplifier 

These  are  linear  push-pull  amplifiers  operated  in  class  B.  (Diagram 
3.2).  The  tuned  circuits  are  heavily  damped  to  obtain  steep  pulse 
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edges.  The  power  amplifier  has  a  maximum  output  of  800  W  unmodulated 
signal  pulses  i.e.  400  W  modulated  signal  pulses.  It  is  gated  at  the 
grids  to  achieve  good  carrier  suppression  between  pulses  and  to  reduce 
anode  dissipation.  In  Figure  3.2  the  shape  and  frequency  spectrum  of 
the  transmitted  pulse  are  depicted. 


prf  s  50e/s 


V»L  ' 


Figure  3 . 2  Transmitted  pulse 


3.8  Pulse  amplifier 


Before  being  applied  to  the  power  amplifier  the 
through  a  pulse  amplifier  (Diagram  3.2).  Here 
and  clamped  by  diodes  to  appropriate  levels  and 
lower  to  the  midtap  of  the  power  amplifier  grid 


gating  pulse  passes 
the  pulse  is  amplified 
fed  via  a  cathode  fol- 
transformer. 


3.9  Power  supply 

Two  stabilized  power  supplies  give  plus  and  minus  300  V  with  very 
little  ripple.  (Diagram  3.3).  The  plus  800  V  source  to  the  power 
valve  screens  is  only  capacitor  filtered  as  is  also  the  4  kV  anode 
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power  supply,  since  these  voltages  only  need  to  be  approximately  con¬ 
stant  during  the  S  msec  pulses  which  are  synchronized  to  the  AC 
mains . 

The  high  tension  supply  is  variable  with  an  auto  transformer.  A  time 
relay  delays  the  high  voltages  allowing  for  a  warm  up  time  of  1  min 
for  the  tubes. 


3.10  Aerials 


Separate  transmitting  and  receiving  vertical  delta  aerials  are  used, 
each  terminated  in  a  1100  ohm  load.  The  two  aerials  are  connected 
to  the  same  28  m  high  vertical  pole,  but  mounted  at  right  angles 
to  each  other  to  minimize  coupling. 


4  DESCRIPTION  OF  THE  RECEIVER  EQUIPMENT 

This  section  describes  the  receiver  stage,  i.e.  the  conversion  of  the 
faint  HF  signal  from  the  antenna  to  two  separate,  amplified  IF  signals. 
The  description  will  be  given  with  reference  to  seven  circuit  diagrams 
(Diagrams  4.1  -  4.7). 


4.1  Antenna  gate 

To  prevent  the  transmitted  pulse  from  reaching  the  receivers  directly 
a  gate  (Diagram  4.1)  is  inserted  between  the  antenna  and  the  recei¬ 
vers.  The  signal  path  is  opened  by  a  negative  pulse  from  the  pulse 
generator.  During  the  negative  pulse  the  two  series  diodes  will  be 
conducting,  whereas  the  shunt  diode  is  closed.  In  the  interval  be¬ 
tween  pulses  only  the  shunt  diode  will  be  conducting,  and  the  signal 
path  is  blocked.  When  the  gating  pulse  coincides  with  the  reflected 
pulse,  only  this  signal  will  reach  the  receivers. 
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4 . 2  Pulse  generator 

The  pulse  generator  (Diagram  4.2)  can  be  triggered  either  from  the 
AC  mains  or  from  an  external  tone  generator  of  p  r  f  in  the  range 
50-200  c/s.  The  monostable  multivibrator  generates  positive  pulses 
which  are  amplified  and  phase  inverted.  The  phase  of  the  output  pulses 
with  respect  to  the  trigger  voltage  can  be  chosen  arbitrarily,  and  by 
means  of  a  potentiometer  the  pulse  width  can  be  varied  from 
0.4  -  4  msec. 

The  negative  pulses  are  fed  to  three  separate  cathode  followers,  one 
supplying  pulses  to  the  antenna  gate  and  the  two  others  supplying  pul¬ 
ses  to  the  IF  gates. 


4.3  LF  generator 

The  LF  generator  (Diagram  4.3)  consists  of  two  oscillators  and  a 
mixer  valve.  The  two  tuning  capacitors  are  mechanically  coupled,  and 
the  output  voltages  from  the  oscillators  are  of  frequency 
f+  =  275  Kc/s  +  f^/2  and  f_  =  275  Kc/s  -  f^/2  respectively.  Care 
has  been  taken  in  design  and  construction  to  avoid  any  interaction  be¬ 
tween  the  two  separate  oscillators.  After  amplification  the  two  sig¬ 
nals  are  combined  in  the  mixer  valve,  and  in  the  output  the  difference 
frequency  f^  appears.  By  a  filter  the  LF  signal  is  extracted  and 
fed  to  a  cathode  follower.  The  output  voltage  can  be  varied  by  means 
of  a  potentiometer. 


4.4  HF  oscillator 


The  frequency  of  the  HF  oscillator  (Diagram  4.4)  can  be  varied  in  the 
range  2.5-6  Mc/s.  The  oscillator  valve  is  followed  by  a  tuned  amp¬ 
lifier  that  also  provides  isolation  between  the  oscillator  and  the  ex¬ 
ternal  load. 
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4 . 5  External  local  oscillator  system 

This  system  (Diagram  4,5)  generates  the  two  mixer  signals  that  sub¬ 
stitute  the  signals  from  the  local  oscillator  in  the  receivers.  The 
two  voltages  are  coherent  and  of  frequencies  f^  =  fQ  +  f^  +  450  kc/s 

and  f_,  =  f  -  f,  ♦  450  kc/s  respectively  where  450  kc/s  is  the 
II  o  L 

frequency  in  the  IF  channel.  The  main  advantage  of  this  system  comes 
from  the  fact  that  the  oscillators  need  not  have  a  high  frequency 
stability.  A  small  relative  shift  in  one  or  more  of  the  generated  fre¬ 
quencies  does  not  affect  the  operation  of  the  complete  system. 

The  frequency  f^  is  generated  by  the  following  process.  The  signal  of 
frequency  f+  is  fed  via  a  cathode  follower  (Diagram  4.3)  to  a  tri- 
ode  pentode  performing  doubling  and  mixing  with  the  hF  signal  of  fre¬ 
quency  fQ  -  100  kc/s.  After  successive  filtering  and  amplification 
the  output  voltage  will  be  of  the  sum  frequency 
fj  =  <fQ  -  100)  +  2 • (275  +  fL/2)  =  f0  +  fL  +  450  kc/s 

Originally  the  frequency  fj  was  generated  according  to  the  equation 
f j  -  fQ  +  2 •  (225  +  f^/2).  This  scheme  failed  because  the  output  from 
the  mixer  in  addition  to  the  desired  signal  contained  weak  signals  of 
frequency  fj.  1  450  kc/s  giving  rise  to  false  signals  from  the  re¬ 
ceiver. 


4.6  Receiver  with  output  cathode  follower 

The  receivers  used  were  the  Eddystone  Communication  Receiver,  Model 
680  X.  Some  minor  modifications  were  necessary  to  incorporate  tiie 
receivers  in  the  complete  system}  (Diagram  4,6). 

1)  The  local  oscillator  valve  was  removed  and  replaced  by  a  connection 
to  the  external  local  oscillator  system. 
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2)  As  the  desired  output  from  the  receiver  is  the  IF  signal,  a  cathode 
follower  was  mounted  into  the  receiver  casing. 

3)  The  receiver  was  provided  with  better  shielding  and  aerial  connec¬ 
tions  to  prevent  any  pick  up  of  the  transmitted  pulse. 


4.7  IF  gate 

To  achieve  extinction  of  the  signal  coming  directly  from  the  transmitter 
it  was  necessary  to  have  an  extra  gate  in  the  signal  path.  This  gate 
was  located  just  after  the  receiver  (Diagram  4.6).  The  operation  is 
almost  identical  to  that  of  the  antenna  gate.  The  resonant  circuit  in 
the  connection  to  the  pulse  generator  provides  a  large  isolation  impe¬ 
dance  to  the  IF  signal,  but  negligible  impedance  to  the  pulses. 


4 . 8  Mains  filter 


A  mains  filter  (Diagram  4,7)  was  used  to  prevent  pick  up  transmitter 
signals  from  reaching  the  receiver.  The  filter  was  inserted  in  the 
AC  mains  line  to  the  receiving  equipment. 


5  DESCRIPTION  OF  THE  DETECTOR  AND  RECORDING  SYSTEM 

This  section  describes  the  processing  of  the  IF  signals  from  the  re¬ 
ceivers.  Although  many  other  more  elaborate  schemes  are  possible,  an 
advantage  of  the  system  used  is  simplicity  and  ease  of  operation. 

The  system  is  described  with  reference  to  five  circuit  diagrams. 
(Diagrams  8.1  -  8.4). 
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.  1  The  phase  changer 

In  order  to  add  the  two  IF  signals  with  different  phase  relationships, 
one  of  the  signals  is  sent  through  the  phase  changer  (Diagram  5.1). 
The  principle  of  this  device  can  best  be  explained  in  connection  with 
Figure  5.1. 


Figure  5.1  Principle  of  phase  changer 

Each  of  the  four  output  voltages  from  the  phase  network  are  of  the 
same  amplitude  and  frequency  as  the  input  ,  but  have  been  shifted 
in  phase  as  indicated.  The  two  voltages  applied  across  the  impedances 
Z  are  thus  in  phase  quadrature.  As  the  impedance  per  unit  length 
is  constant,  the  signal  picked  up  by  the  sliding  contact  is  propor¬ 
tional  to  the  distance  a  from  the  mid  point. 
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The  movements  of  the  contacts  are  coupled  in  such  a  way  that  the  two 
deviations  a,  and  a2  are  also  in  phase  quadrature  (see  Figure 
5.1  b).  The  output  voltage  VQ  will  then  be  identical  to  the  input 
voltage  except  for  a  phase  difference  8  that  is  determined  by  the 
positions  of  the  sliding  contacts.  When  the  mechanical  movement  cor¬ 
responds  to  a  uniform  rotation  f  in  Figure  5.1  b  ,  the  phase  8 
will  be  a  linear  function  of  time.  The  effect  on  the  output  signal  is 
thus  equivalent  to  shifting  the  frequency  slightly  by  an  amount 

f  =  de/dt. 

r 

The  details  of  the  phase  changer  are  given  in  Diagram  5.1.  The  input 
signal  is  amplified  and  transformer  coupled  to  a  phase  network  consist¬ 
ing  entirely  of  resistors  and  capacitors.  The  operation  of  this  phase 
network  is  independent  of  frequency  in  the  vicinity  of  f .  .  By  means 
of  potentiometers  the  amplitudes  and  phases  of  the  output  signals  can 
be  given  the  proper  values. 

The  impedance  Z  is  a  capacitor  consisting  of  an  inner  cylinder  and 
an  outer  cylinder  divided  into  two  parts  electrically  insulated  from 
each  other.  The  inner  cylinder  is  free  to  move  and  is  equivalent  to 
the  pick  up  contact  in  Figure  5.1.  The  two  inner  cylinders  are  coupled 
and  driven  back  and  forth  by  an  electro  motor,  the  device  very  much 
resembling  a  steam  engine. 

The  signal  currents  from  the  pick  up  cylinders  are  added  and  amplified 
in  a  tuned  amplifier  with  variable  gain. 


.  2  Addition  circuit  and  linear  detector 


The  two  echo  signals  are  combined  and  detected  in  the  circuit  shown  in 
Diagram  5,2.  The  summation  is  performed  by  an  amplifier  with  large 
negative  feedback.  To  reproduce  the  envelope  of  the  pulses  properly, 
the  time  constant  t  of  the  detector  must  satisfy:  l/f^<<  t  <<  T  , 
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where  T  is  the  pulse  width.  Due  to  the  phase  changer  the  amplitude 
of  the  video  output  pulses  will  be  oscillating  with  period  fr. 


5 . 3  Box  car  detector 

The  box  car  detector  (Diagram  5.3)  is  a  peak  voltmeter  with  very 
small  time  constant  for  both  charging  and  discharging  when  working  on 
pulse  signals.  On  the  other  hand,  the  time  constant  in  the  interval 
between  pulses  is  approximately  infinite.  This  behaviour  of  the  peak 
voltmeter  is  very  important  when  the  input  signal  is  varying  rapidly. 

The  desired  characteristics  are  obtained  by  replacing  the  usual  dis¬ 
charging  resistor  in  a  normal  single  diode  rectifier  with  a  valve. 
This  is  biased  beyond  cutoff  between  pulses,  while  it  is  fully  con¬ 
ducting  when  the  signal  pulses  are  arriving.  In  Figure  5.2  the  wave¬ 
forms  in  the  box  car  detector  are  shown. 


Input 

pulses 

prf  =  50c/s 


Trigger 
pulses 
pr »  *  50c/s 


Output 
vol  tage 


Figure  5.2  Principle  of  box  car  detector 
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To  achieve  proper  operation  of  the  box  car  detector,  the  trigger  pulses 
should  have  a  minimum  height  and  be  of  somewhat  shorter  duration  than 
the  echo  pulses  as  indicated  in  the  figure.  The  output  signal  is  then 
the  envelope  of  the  input  superimposed  on  a  DC-voltage, 


5 . 9  Trigger  pulse  circuit 


The  trigger  pulses  are  generated  from  the  echo  signals  in  the  nonlinear 
circuit  shown  in  Diagram  5.4.  By  successive  amplification,  detection, 
clipping  and  shaping,  the  output  pulses  are  rendered  almost  independent 
of  the  characteristics  of  the  echo  pulses. 

The  signals  in  the  two  IF  channels  are  amplified  and  detected  separate¬ 
ly.  After  addition  the  video  pulses  are  again  amplified  and  then  fin¬ 
ally  clipped  and  shaped.  Depending  on  the  amplitude  of  the  echo  pul¬ 
ses,  the  valves  in  the  circuit  will  operate  as  linear  amplifiers 
(small  signals)  or  limiters  (large  signals). 

The  magnitude  of  the  pulses  are  governed  by  the  two  resistors  in  the 
voltage  divider.  The  width  and  shape  of  the  trigger  pulse  are  to  a 
certain  extent  determined  by  the  components  in  the  pulse  shaper. 


5 . S  Power  supply 

The  power  supply  serving  the  detector  and  recording  system  is  shown  in 
Diagram  5.5.  It  provides  a  stabilized  voltage  at  300  V  plus. 


5.6  Brush  recording  equipment 

The  recorder  used  is  a  two  channel  Brush  recorder  Type  BL202  equipped 
with  a  DC  preamplifier.  The  speed  of  the  paper  chart  is  normally 
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5  mm/s  ,  but  in  some  cases  a  higher  speed  has  been  used.  The  recording 
on  paper  chart  is  used  mainly  for  adjusting  and  supervisory  purpose. 
Figure  3.4,  Part  II  shows  a  sample  of  a  typical  recording. 

The  signal  is  recorded  in  the  lower  channel,  whereas  each  mark  in  the  . 
upper  channel  corresponds  to  a  sample  taken  by  the  A-D  converter. 


5 . 7  Analog-digital  converter  and  tape  recorder 

It  is  possible  to  extract  the  desired  information  from  the  recording 
on  paper  chart,  and  in  the  first  phase  of  the  experiment  this  was 
actually  done.  But  the  manual  analysis  of  even  a  few  minutes  regi¬ 
stration  is  very  tedious,  and  facilities  for  recording  on  digital 
form  were  therefore  added  to  the  equipment.  The  A-D  converter  and 
tape  recorder  are  parts  of  a  data  processing  system  developed  at 
HDRE . 


The  sampling  times  are  controlled  by  the  electro  motor  that  drives 
the  phase  changer.  As  is  seen  from  Figure  3.4,  Part  II,  four  samples 
are  taken  per  revolution  of  the  phase  changer.  Four  protruding  knobs 
spaced  90°  apart  on  a  wheel  trigger  the  sampling  process.  Each 
sampled  value  is  converted  to  a  nine  bit  number  and  finally  recorded 
on  magnetic  tape. 


6  ADJUSTING  OF  EQUIPMENT 

when  measurements  are  made  some  care  must  be  exercised  to  obtain  pro¬ 
per  operation  of  the  equipment  and  these  will  be  mentioned  briefly. 

It  is  desirable  to  have  a  transmitted  pulse  of  long  duration,  and  the 
pulse  width  should  consequently  be  close  to  the  maximum  set  by  the 
height  of  reflection.  Furthermore  the  balanced  modulator  in  the  trans¬ 
mitter  should  be  adjusted  to  get  complete  suppression  of  the  carrier. 
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The  simultaneous  tuning  of  the  receivers  is  straightforward,  but  care 
must  be  taken  to  provide  equal  gain  in  the  two  IF  channels.  This  is 
■easy  to  adjust  when  the  transmitted  pulses  are  received  directly  (with 
sufficient  attenuation).  Due  to  the  phase  changer  the  recorded  signal 
will  be  oscillating,  the  minimum  being  zero  when  the  gains  are  equal. 

It  is  essential  that  only  the  wanted  pulses  pass  through  the  gates,  as 
any  other  signals  will  currupt  the  operation  of  the  box  car  detector. 

The  position  and  width  of  the  gating  pulses  must  therefore  be  adjusted 
carefully. 

After  adjusting  the  equipment  according  to  these  rules  the  measure¬ 
ments  proper  are  very  easy  to  perform,  and  there  is  no  need  for  frequent 
readjustments . 
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Diagram  4 . 7  Mains  filter 


Diagram  5 . 1  Phasa  changer 


AMPLIFIER  PHASE  NETWORK  AMPLIFIER 


ADDING  CIRCUIT  CATHODE  FOLLOWER  LINEAR  DETECTOR 


TRIGGER  PULSE  IN 


I 

CO 

CV*  o 
10  M 
O  CO 

H  £ 
to  a, 
o 

fu  05  w 
<  fn  cj 


OS 

5 

X 

o 

u 

OS 

OS 


G  X 

kO 

3- 

a 

w 

os 

u 

£ 

CO 

2 

o 

u 

03  < 

Q 

2 

CO 

G  <  3 

2  rrj 

o 

OS 

O  H  as 

to 

M 

4 

2  CO  O 

H  cn 

X 

MUZ 

OS  P 

a 

&4 

Pi 

O 

o 

O  CO 

y  x  - 

CU  G 

H  as 

Q  o  as 

C4  a 

X 

co  o 

as  t4 

os  p 

a 

HPZ 

_  <  G 

0 

3 

os  u  u 

3  (4  G 

14  £ 

H 

X  <TSl 

a:  co  w 

OS 

CO 

u  re  m 

o  u  ^ 
x  as  !*S 

Q  O 

2  P 

< 

Ou  E-*  X 

0 

• 

o 

C  3 

0  0 

o 

O  0  O 

G  >v 

to 

<H  P  *H 

H  CO  0 

1 

P  0  3 

P  X 

»-> 

CO 

O  000 

CN 

CJ 

e-H  p 

•  CO  o 

0 

M 

p 

CO  C  o 

o 

CO 

P  0  c 

0-P  H 

\_s 

X 

©  0  0 

•H  ITJ 

P 

X 

P  >  0 

og 

to 

a 

3  C  5 

c  a)  3 

o 

OOP  P 

a)  a  c 

Cm  OS  (4 

0  0 
&g* 

3  0 

0  o 

<  H  O 

•H  fl  C  P-H  O 
o  *PP  H 
3  ®P  10 
G  P  P  3P 
0  0  0  0  V  C 
GP  OP,  0 
0  a  0  0  u  0 
P  0  P  OO  2 
3  0  P  W  U  P 
P  C  0  0 

o  o  O  >X  A 
3  P  HO 
P  0  ©>H  0 
P  4I£  U)X  0 
0  G  P  O  £  p 
P  P  P 
P  P  OU 
0  «P  0  0  3 

P  0  P  >  C 
3  >>  0  O  fl 

00  3 

0  P  3  ■O  0  0 

u  0  p  a>  aop 

P  X0  PGA 

-4  0  CJ  0  0 

H  II  0  U  p 

H 

P  >♦  0 
MU  > 
P  G 
0  X 
co  3  P 
0  CTA 
>  0  OO 
0  P  P 
3  <M  A 


i  g  p 


H 

O  G 

2  0  OO 
G  3 
H  P  O 
CO  P 
00  C  G 
<  P  P 


G  X 
CO  < 
u  <36 
O  H  as 

2  CO  O 

uuz 

p 

y  x  * 
Q  CJ  as 
OS  (4 
<  G 
S  W  G 
as  co  u 
o  w 


co 

C4 


as 

3 

S3 

CD 

C4 

as 

as 


as 

(4 


o 

2  ro 

to 
H  cn 

as  p 
o 

04  G 
(4  O 
as  p 
0 

W  £ 

as 

a  o 

2  p 


CO  • 

G  0  (0 

O  3  0  G  ^ 

P  0  P  H  0  O 

PPT3  P  M 


O  0  0 

e  oo  p 

•H 

PPG 
0  CO  0 
P  0  0 
3  >  J 
OO  C  P 
0 
p 


•  wo 
CO  c  o 
0., 

H  i 

o  u 

•3 

c 

0 


sr 


H  0  3 

G  a4  c 

2  0  0 
C  P  P  o 
0  O  <P  P  p 
•H  0 
0  P  3  P 
P  0  0  0 

OP 


0  G  0  0  p.  i 
pap-* 

3  w  p 

poo 


2 

14 

CO 

as 

M 

X 

a 

o  co 
H  as 
co  o 
m  a*  2 
as  o  w 

X  <SL 
Olffl 

a  h  x 


c 
0 

-  ..I  0 

a  o  3t 

co  a  p 

0 

O  XG  G 

3  0  r-4  o 

p  P  0  0  p  CO 

P  G  CO  G  0 

W  0  P  0  3  P 

GPP 
P  P  U-4  O  P 
0  O  0-0 
P  <P  P  >  C 

3  0  X  0  O  0 

oo  3 

2  P  3  3  0  0 
0  P  0  OOP 
P  X  CO  P  G  G 

W  0  O  0  0 

P  0  0  Pp 
I  GP  P 
..  U,  p  in  X  0 

H  0  O  > 

O  G  P  C 

2  0  OO  0  >> 

G  3  CO  3  p 

h  t*  o  «  cr.c 

«  >i  >  V  M 
M  C  X  4  C  -H 
<  -H  « 


m  p 

O  CO 

^  X 

P  3S 

CD  (L, 

o 

a  as  u 

<H  c 5 


G  X 
CO  < 
W<3 

oho: 

2  CO  O 
MUZ 

U4 

Ml  * 

Ou  0! 

3a 

3  14  J 
«COW 
014  *-a 

20S  *S 


OS 

u 

X 

a 


W 

• 

o 

C  3 

0 

CA 

o 

O  0  0 

G 

V 

to 

•H  PP  H  0 

0 

1 

P  0  3 
0  OO  0 

P 

0< 

CN 

CO 

CJ 

S  P  P 

•  0 

o 

to 

p 

P 

w  e 

o 

O 

CO 

P  0  C 

0  p 

p 

w 

X 

0  0  0 

2 

P 

X 

P  >  0 

y  a 

to 

a. 

3  C  5 

C  V 

3 

o 

OOP  P 

2  u 

C 

u. 

as  (4 

0  0 
P  0  G 
P  P 

P  0  C 

3 

cr  c 
0  o 

P  p 

0 

o 

<  H  O 

M 

0  4 

p 

H 

3 

0P  _  0 

M 

c 

P  P  3  H 

as 

0 

0  0  0  9 

c 

GP  0  h 

0 

G 

0 

a  0  0  u 

0 

X 

p 

w  p  a  o 

3 

CJ 

3 

o  p  «  o 

p 

X 

P 

c  0 

0 

as 

O 

0  O  XA 

G 

as 

3 

P  p  0 

r,  w  -  H  „ 

P  0  G  W  A  0 
0  G  P  0  J  H 


P  «H  O 
0  *P  O 
P  O  P 
3  X  0 


4- 

CO 

2 

00  3 

o 

W 

©  p  3  3 

0 

2 

0 

P  0  P  0 

2  ro 

o 

G 

P  X  0  P 

to 

M 

U 

M  0  0 

H  cn 

X 

P  0  0 

06  P 

(M 

Ck 

1  G  H 

O 

O 

O  CO 

..  a  p  »h 

a  g 

h  as 

H  0 

(4  a 

X 

C5  O 

CJ  G  fc 

as  p 

a 

H  P  2 

<  0  DC 

0 

X 

G  CJ  M 

OS  G  3  0 

W  £ 

Eh 

x  <  ta 

H  P  0  0 

OS 

Q  O 

CO 

0  E  O 

CO  P  > 
CO  CG  0 

2  P 

< 

a  H  x 

<  P  P  3  ' 

h 

0  3 
>  C 
0  0 

0  0 
OOP 


0  > 

c 

2  * 

3  P 


G>« 
m< 
M  <3 
o  h  as 
2  CO  O 
U  14  2 
P* 

MI  • 
o  oos 


3  14  J 
OS  CO  (4 
O 

2  as  *s 


O 

2  <*> 
to 
H  o> 
as  p 

o 

a.  g 

C4  O 

os  p 

a5 

Q  o 

2  H 


as 

3 

X 

CD 

(4 

OS 

as 


os 

(4 

£ 

CO 

o 

2 

o 


2 

U 

CO 

OS 

(4 

£ 

H  OS 
COO 
HPZ 

w 

Olo 

a  h  x 


CO  • 

c  3  0  0 

0  0  O  G 

p  p  p  H  a  u 

a  &*2  ’H* 

0  00  0 

S  *H  P  •  0  O 
,  +■  a  c  o 
P  0  C  0-H  P 
0  0  9  -H  0 
p  >  «  a  a 
3  G  3  **  “  - 
00‘H  P 
0 

Eg- 

•H  0  c  £  O 

o  <H  P  P 
■O  0 

GPP  3  P 

0  0  0  0  0  c 

G  P  O  p  0 

v a®  0  p 0 
p  0  P  ao  i 
3  O  P  «  o  P 
P  c  o  0 
0  o  o  >wcZ 

3-H  P  cJ 

P  0  ©*H  0 
P  ©  G  0  G  0 
•  X  P  0  J'H 


c  0  T3 

•  *  C  c 

:  g.® 


P 

0  <P  O 
H  O 
3  >> 

00  TJ 

9  P  3 
P  0  P 
P  >t0 

H  0 

P  0 


•mop 


0  O  0 

T3  0  0 

Sg-3 

y  -g?* 

<  0  00 

OS  X  3 

h  p  r 

co 


g  cx 


0  >s 
0  3  P 
0  0*G 
>#M 
0  P  -H 
G 


